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1. Introduction 

1.1. Motivation 

Unmanned Underwater Vehicles (UUV’s) perform many of their missions in shallow 
water environments subject to the forces of ocean waves and the proximity to the ocean 
floor. Under these conditions, accurate vertical position control is necessary to prevent 
broaching or hitting the ocean floor. Accurate horizontal position control is necessary to 
enable the UUV to conduct its mission with accuracy and return to a predetermined 
recovery point. Shallow water position control is made more difficult by ocean waves. 

In deep water the effects of these waves are negligible, but the effects in shallow water 
are significant. Important shallow water missions include pollution monitoring, marine 
life sampling, bottom contour mapping, and mine location. 

Currently, UUV’s are controlled in shallow water by altering empirical control 
parameters for better shallow water performance and by establishing empirically based 
operating depth limits on the UUV operations. These operating depth limits are based 
upon wave conditions. With a thorough understanding of the dynamics of UUV’s in 
shallow water and the forces and moments on vehicles due to sea waves in these waters, 
improved control systems and vehicle designs can be achieved to allow the UUV to 
operate in shallower water and in larger waves than is commonly done. This will allow 
the UUV to be more effectively perform its missions. 

This thesis explores the effects of variation in water depth and vehicle submergence 
on added mass, damping, and restoring forces. 

1.2. Historical Background 

This work builds on the work that Timothy Prestero performed to build a 
mathematical simulation of the REMUS behavior in deep water. This work is reported in 
his Master of Science Dissertation for the Massachusetts Institute of TechnologyAVoods 
Hole Oceanographic Institution Joint Program in Oceanography/Applied Ocean Science 
and Engineering entitled “Verification of a Six-Degree of Freedom Simulation Model for 
the REMUS Autonomous Underwater Vehicle.” Mr. Prestero calculated the 
hydrodynamic and hydrostatic coefficients based upon deep water performance far from 
a boundary surface. This thesis extends Mr. Prestero’s work by determining those 
coefficients in shallow water and near the surface. ’ 

1.3. Research Platform 

The platform for this research is the REMUS (Remote Environmental Monitoring 
Units) AUV (autonomous underwater vehicle) developed by the Oceanographic Systems 
Laboratory at the Woods Hole Oceanographic Institution. This low-cost, modular AUV 
was developed for coastal monitoring and multiple vehicle survey operations.^ REMUS 
has also been adopted for use in mine-counter measure operations for the United States 
Navy.^ REMUS has most recently been used by the United States NaAry to hunt for 
mines from the Iraqi port of Umm Qasr in support of Operation Enduring Freedom.'* 




1.4. Assumptions 

To simplify analysis, the author made the following assumptions: 

• The vehicle is port-starboard symmetric. 

• The vehicle is a rigid body of constant mass. 

• There are no significant vehicle dynamics occurring faster than the data sampling 
frequency of 25 Hz. 

2. The Coordinate System 

The coordinate system used for this research is shown in Figure 1. This is a body- 
fixed right-handed coordinate system with the x axis defined along the axial length of the 
vessel and the z axis defined downward. The origin of the body-fixed coordinate system 
is at the vessel amidships. The variables shown in Table 1 are defined using the 
coordinate system shown in Figure 1. 




Figure 1. Sketch showing positive directions of axes, angles, velocities, forces and moments. 
(Feldman, 1979) 


Table 1. Coordinate System Variables. 


X 

y 

z 

u 

V 

w 

p 

q 


Surge position forward. 

Sway position to the right. 
Heave position downwards. 
Velocity in the surge direction. 
Velocity in the sway direction. 
Velocity in the heave direction. 
Rotation about the x axis. 
Rotation about the y axis. 




r Rotation about the z axis. 
X Force in the x direction. 
Y Force in the y direction. 
Z Force in the z direction. 
K Moment about the x axis. 
M Moment about the y axis. 
N Moment about the z axis. 


3. The Equations of Motion 

This work primarily explores the forces and moments in sway, heave, pitch and yaw 
due to motion of the vehicle itself in finite depth water. For each series of tests, the 
vessel was moved in only one plane at a time. The resulting hydrodynamic forces were 
determined by subtracting inertial forces from the measured forces. Then, the 
hydrodynamic coefficients were extracted fi'om the hydrodynamic forces. The equations 
of motion were used to perform the mathematical operations. 

The forces and moments associated with surge and roll have not been investigated. 

3.1. Vessel Inertial Dynamics 

The linearized equations of motion with a body-fixed coordinate system for an 


unrestrained vessel in water are given by 

X = m\ 

[ii-vr + wq- x^q^ + r^) + y^ipq - r) + z^ipr + q)j 

Y = m\ 

[v-wp + ur-y^^ir^ + p^) + Z(^iqr-p) + Xc(qp + r)] 

Z = w| 

[w-uq + vp- Zcip^ + r^) + x^(rp -q) + y^irp + p)] 


K = 1^P + (4 - Iyy)^r - I^ipq -I- r) + 4(r^ - q^) + I^ipr -q) + 

m[y(;(w + pv-qu)-Zc(v + ni-pw)] (1) 

M = I^q + (4 - 4)/7r - 4(p + qr) + IJp^ - r^) + I^(qp -r) + 
m[Zg{u + qw — rv) — XQ(w + pv-qu)] 

N = IJ + (l^-IJpq-I^iq + rp) + l^^Xq^-p^) + I^^(rq-p) + 
m[X(;(y + ru —pw) —y^(u + qw-rv)] 

where 

m is the mass of the vessel 

(xG,yG,ZG) are the coordinates of the center of gravity of the vessel in the body 
fixed coordinate system. 

Ijk are the moments of inertia. 

These equations can be simplified by fixing the coordinate system at the midship 
location of the vehicle. The equations can also be simplified by assuming that the lateral 
distance from the midship location to the center of gravity is negligible, i.e. yG =0. 

Further simplification can be obtained by testing and analyzing motions in the vertical 




and horizontal planes separately. This research does not examine hydrodynamic forces in 
surge and roll, so the relevant simplified equations are: 

Y = m [v + t/r + x^r] 

Z = m^w-Uq-Zfjq 

M = Iyj,q + m [zg {u-vr + wq)-Xf^{w- Uq)] 

N = IJ- + tyiXq (v + rw) 



3.2. Hydrodynamic and Hydrostatic Equations 


This thesis explores the hydrodynamic forces and moments due to unsteady motion of 
an underwater vehicle. For that reason, hydrostatic effects have been removed from the 
data by subtracting the mean forces and moments from all measured forces during the 
analysis. 

The forces and moments experienced by a ship are assumed to be the forces and 
moments arising from motions of the ship which in turn have been excited by another 
source. These forces and moments are computed as functions of speed and acceleration. 
A mathematically useful form is derived using the Taylor expansion of a function of 
multiple variables. For example, sway force, Y, and yaw moment, N, are represented 
functionally as 


Y = Fy{u,v,u,v,r,r) 

N = FXu,v,u,v,r,r) 

The Taylor expansion of a single variable states that if the function of a variable, x, 
and all its derivatives are continuous at a particular value xi, then the value of the 
function at a value of x close to xi can be expressed as 






dx 


dx^ 




. dVix) 


dx^ 


dx” 


(4) 


where * 


f(x) is the value of the function at x close to xj 
f(xi) is the value of the function at x = Xi 

(^X = X - Xi 


d"f{x) , 


dx” 


is the nth derivative of the function evaluated at x = X] 


By making 6x sufficiently small, higher order terms can be neglected. Equation (4) 
reduces to 

df(x) 


f(x) = fixi) + Sx- 


dx 


(5) 


and is called the linearized form of the Taylor expansion. 

For functions of two variables the linearized form of the Taylor expansion is 

OX dy 

Again, (5x and ^y must both be small enough that higher order terms can be neglected. 


(6) 




Hydrostatic motion stability typica;lly considers the effect of very small perturbations 
on the behavior of the ship. Thus, the linearizing assumption for the Taylor expansion 
can be used to describe the hydrodjmamic behavior of a body. Analysis of data from this 
research indicates that similar non-dimensional results were obtained for tests done at 
different amplitudes and the Fourier coefficients at the excitation frequencies dominated 
all others. Because of these facts, the linear terms do indeed predominate and the model 
based on them is sufficient to describe the relation between vehicle motions and the 
forces and moments they generate. Using the linearized Taylor expansion, equation (3) 
can be written as 

Y = F(u^,v„u^,v^,^,r^) + iu-u^)— + {v-v^)—+... + (r-r^)-— 

du ov dr 

A = F,(Mi , V„M,, v„ r,, + (m - M,)—+ (v - V,)-—+...-I-(r - r, )— 

du dv dr 

At this point several simplifying assumptions can be made. The first assumption is 
that the initial motion is in a straight line at some constant speed. Therefore, 
ti, =Vi=^;=7^=0. The ship is symmetrical about the xz-plane, so v/ = 0. Symmetry 

also leads to the conclusion that dYldu = dY ldu = 0 because forward motion will not 
cause a lateral velocity. Also, a ship traveling forward in equilibrium in straight line 
motion experiences no sway force, so the term Fy{u^,v^,u^,v^,r^,r\) is also zero. The term 

U] is equal to the straight line velocity U. These assumptions reduce equation (7) to 

^ dY dY . dY dY . 

dv dv dr dr 

,, dN dN . dN dN . ^ ^ 

dv dv dr dr 

In the simplified notation used by the Society of Naval Architects and Marine 
Engineers and including Pitch and Heave, the simplified linear hydrodynamic equations 
become^ 

y = V+y,v+y^r+y/ 

N = Nv 4 - N-v + Nr + N-r 

(9) 

Z = Z^w + Z^w+Z^g + Z^q 

M = M^w + M^w + M^q + M^q 

The simplified notation is interpreted such that Y^v is the sway force related to sway 
motion and Y^ is the maneuvering coefficient of sway force due to sway motion. 

In accordance with the standard notation the terms of equation (9) include the effect of 
the rudder and stem planes held at zero degrees. The experiments to extract the 
coefficients were all performed with no deflection of the control surfaces. Other 
experiments were performed with control surface deflection. For those experiments, 
equation (9) has additional terms related to rudder and stem plane angle.^ 

3.3. Added Mass and Damping 

The hydrodynamic forces relating to the motion of the body in the fluid can be divided 
into components in phase with the acceleration and components in phase with the 
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velocity of the body. The hydrodynamic force due to the acceleration of the body in a 
fluid is known as an added mass force. The hydrodynamic force due to the velocity of 
the body in the fluid is known as a damping force. These forces are can be discerned by 
their phases relative to the driving motion. Forces in phase, but opposite in sign, with the 
driving motion are related to acceleration and are added mass forces. Forces 90 degrees 
out of phase with the driving motion are related to velocity and are damping forces. In 
terms of complex notation, the added mass is related to the real component of the 
measured force and the damping is related to the imaginary component of the measured 
force. 

4. Non-Dimensionalizing 

Throughout this thesis several quantities are given in both dimensional and non- 
dimensional form. Final results are given in non-dimensional form to be readily available 
for use with other bodies of similar shape. Non-dimensional quantities are denoted by a 

7 

prime symbol ('). The equations for non-dimensionalizing are: 

Y' = - - - 




\pW 


\pU^L^ 


\pW 
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M 


M 


\pU^L^ 


1^4 ' ^ 

M... = 


\pUL^ 


A4 ' ^ 

M.:. = 


hpl" 


M’ =-^ 

’ i pW 

, M, 

M =—V 
’ \pL^ 


N 


N 


K. 


N, = 


\pu^n 

K 

\pVl} 

K 

kpL' 

\pW 

N'=-^ 

' \PL' 


Other non-dimensional equations include 



m - 


m 


\pL^ 


/ ’ = 




\pL^ 


~ L 
U 

, 1 

V - —V 

U 

., L . 

V = — 


w — —w 
U 

., I . 

w = —-w 
L 

q = —q 
U 

■ H • 




L 

r = —r 
U 

. 1 } . 

r = —-r 


Fr = 


CO = 


U 


CO 


Submergence' = • 


Length 


Submergence 

Using non-dimensional coefficient, the equations of motion have the form 

y = 7 'v'+y 'v'+7 'r'+Y 'r' 


( 12 ) 

(13) 


5. Experimental Procedure 

The determination of the maneuvering coefficients was conducted using both full 
scale and model scale experiments. Full scale experiments were used to determine the 
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body lift and control surface effects. Model scale experiments were used to determine 
the unsteady motion effects. 

5.1. Experiment Apparatus 

5.1.1. Model Geometry 

Figure 2 shows the geometry of the full scale model. The small scale model is 
geometrically similar at a scale of 0.4334. This scale was selected to provide the smallest 
model that would contain the transducer discussed in Section 5.1.2 without incidental 
contact between the transducer and the model. 



Ballast 


Forward 

Transducer 


Figure 2. Full Scale Model Geometry 

Figure 3 shows the full scale model mounted in the United States Naval Academy 
Towing Tank. Figure 4 shows the 0.4334 scale model mounted in the Massachusetts 
Institute of Technology Marine Computation and Instrumentation Laboratory. 







Figure 3. Full Scale Model Mounted in United States Naval Academy Towing Tank 



Figure 4.0.4334 Scale Model Mounted at the MIT Marine Instrumentation and Computation 
Laboratory 

5.1.2. Force and Moment Measurement 

The forces and moments were measured using a UDW3 underwater transducer 
manufactured by Advanced Mechanical Technology, Inc. The transducer, shown in 
Figure 5, is able to simultaneously measure forces and moments in all of the three 
orthogonal directions (making six measurements of forces and moments) and is suitable 
for underwater applications. A pressure compensating bladder in the transducer equalizes 
internal and external pressures to allow underwater operation with little effect of 
hydrostatic pressure. The capacities and general specifications of the dynamometer are 
shown in Table 2. 
















Figure 5. UDW3 Underwater Sensor.* 


The transducer was mounted to a bulkhead within the volume of the vehicle. The strut 
was attached to the end of the transducer not attached to the vehicle. Sufficient clearance 
was provided to ensure the transducer output was not compromised by contact with the 
sides of the vehicle. 


Table 2. Dynamometer Capacity and Specifications’ 


Vertical and Lateral Force Capacity 

556 N 

Axial Force Capacity 

1112.1 N 

Pitch and Yaw Moment Capacity 

28.2 N-m 

mmssmmmmmm 

14.1 N-m 

Vertical and Lateral Force Sensitivity 

2.7 /uVI{V*N) 

Axial Force Sensitivity 

.61 iuVI{V*N) 

Pitch and Yaw Moment Sensitivity 

m.2^iV l{V*N-m) 

Roll Moment Sensitivity 

97.4/iF/(F*V-m) 

Vertical and Lateral Force Stiffness 

5.3x lO^’N/m 

Axial Force Stiffhess 

7.88 X 10' N/m 

Roll Moment Stiffhess 

5.7 X 10^ N-m/radian 

Weight 

2 kg 

Recommended Excitation 

10 V or less 

Crosstalk 

< 2% on all channels 

Temperature Range 

-17 to 52° C 

Force Channel Hysteresis 

± 0.2% Full Scale Output 

Force Channel Non-Linearity 

± 0.2% Full Scale Output 


Excitation for the transducer and amplification for the output were provided by a 
MSA-6 Mini-Amplifier also developed by AMTI. This amplifier, shown in Figure 6, 
provides excitation and amplification for up to six channels. The excitation is selected by 
individual jumpers for each channel and ranges from 2.5 to 10 volts. The gain for each 
channel is also selectable by jumpers and ranges fi'om 1000 to 4000. The output of the 
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amplifier is ±10 VDC. The amplifier contains an auto-zero feature that allows for push 
button zeroing of the output of the load cell. 



Figure 6. MSA-6 Mini Amplifier’* 

The output of the amplifier was connected to an analog-to-digital converter installed 
in a notebook computer. The system control software sampled the six channels of output 
of the load cell and the six positions of the gantry system at an operator selected 
frequency of 25 Hz. 

5.1.3. United States Naval Academy Tests 

Full scale model testing was performed at the United States Naval Academy 
Hydromechanics Laboratory shown in Figure 7. This set of tests included determining 
the forces and moments resulting from body angles in pitch and yaw and control surface 
angles. The towing tank used was 120 ft long, 8 ft wide, and 5 ft deep. The towing tank 
included a wave making machine, a wave absorbing beach and a moving carriage.'^ 



Figure 7. United States Naval Academy Hydromechanics Laboratory Towing Tank 





# 



m 




5.1.4. Massachusetts Institute of Technology Tests 

Small scale model testing was performed at the Marine Instrumentation and 
Computation Laboratory at the Massachusetts Institute of Technology. These 
experiments were performed in order to determine the forces and moments associated 
with unsteady motion. The model was moved in prescribed sinusoidal motions and the 
resultant forces and moments were measured. 

The laboratory contains a tank and a gantry system capable of simultaneous motion in 
five degrees of freedom. The experimental tank is 10 m long, 4 m wide, and 1 m deep. 
The gantry system consists of five different motors and several gear assemblies to ensure 
smooth operation at the speeds and frequencies required for the experiments. The gantry 
system is computer controlled for precise positioning. 

The testing and correction of the very sophisticated gantry system and control 
software occupied a significant amount of the time allocated for the performance of this 
research. The gantry and control software was designed and assembled by D’Ambra 
Technologies, the only firm known to the research supervisor to be capable of developing 
the system and the software. The original contract called for completion of the gantry 
system by January 2002. 

Testing began in June 2002 and significant problems with the system and Control 
software were soon identified. Correction of the problems introduced several weeks of 
delay. The cycle of problem identification and correction continued until very early in 
2003. In this process the vertical axis controls were completely redesigned. The original 
stepper motors were found to be inadequate and were replaced by servo motors. The 
pitch mechanism was strengthened three times to be able to provide the desired frequency 
and amplitudes of oscillation. The gantry system and control software were believed to be 
reliable and accurate in early April 2003. 

The research team also encountered problems related to unidentified faults in the force 
measurement system. AMTI conducted significant troubleshooting of the load cell and 
connections on several occasions to determine the cause of the abnormal readings. Some 
of the abnormal readings were attributed to a fault in the cabling and others were 
attributed to the high level of electronic noise in the long cables of the system. 

Eventually, all of the issues with the force measurement system were corrected. 

One difficult issue with the force measurement system that was discovered late in the 
process is that the measured forces and moments often represent a very small fraction of 
the capability of the transducer. The manufacturer states that the transducer provides 
accurate results at very small fractions of its capability, but this needs experimental 
verification. 

Several hundred experiments were conducted during the process of identifying and 
correcting system problems. Experiments were performed by the author, by an 
independent contractor, and by several undergraduate students acting with supervision. 
The data from these experiments needs to be closely examined to determine if the 
experiments are valid. The first future work that will be done is to predict what those 
tests should show and check for agreement. If the data are found to be valid, they will 
contribute to a more complete data set with less variance that will allow for better 
modeling of the vehicle behavior. 
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5.2. Design of Experiments 

The variables that affect the behavior of an underwater vessel include the depth of the 
water, the submergence of the vessel, the forward velocity, and the frequency and 
amplitude of oscillation. The Central Composite Method was used for Design of 
Experiments in order to reduce the total number of experiments required. 

The Central Composite, or Box-Wilson, Design is a three- or five-level design that 
includes the comer, center, and axial points of the design space. The three-factor Central 
Composite Design space is shown in Figure 8. 


Velocity 



Figure 8. Central Composite Method 


The three factor design space is developed from 15 point designs; a center point 
design, eight comer point designs, and 6 axial point designs. This model represents the 
response surface more accurately than most other methods since the comer points are 
included. Comer points represent the limits of the experimental space. However, 
attempting to reach these comer point designs may strain the engineering model’ 

The selection of test points in an incomplete matrix on the basis of orthogonal numeric 
functions is fine when the dependent variable depends linearly on the input variables. 
However, for things like a nonlinear relation between force coefficient and excitation 
frequency, it is better to be sure that all comers in the test space are tested so that the 
mathematical model will interpolate rather than extrapolate. 

Several experiments that had been planned were not performed due to limitation of the 
gantry system at higher speeds and higher frequencies of oscillation. Other experiments 
were not performed due to physical constraints of the gantry system. Appendix A lists 
the full scale experiments that were performed. Appendix B lists the model scale 
experiments that were used for analysis. All submergences listed in the test matrices are 
to the center of the body. 

A large number of other experiments were performed earlier in the test program, but 
uncertainty in the equipment behavior resulted in uncertainty in the quality of the data 
and data from those experiments was not used. That data will be reevaluated as part of 
future work. 

5.3. Analysis of Experimental Data to Extract Measured Forces and 
Moments 

For testing performed at MIT, the conversion from raw force, moment, and position 
data was performed using MATLAB routines developed by the author. The transducer 
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provided data on forces and moments in the form of voltages for each channel that had to 
be converted to the MKS system. The gantry system provided data on the position of the 
system in a numerical format that had to be converted to the MKS system for analysis. 

For testing performed at the United States Naval Academy, raw force and moment 
data were converted using routines developed the author and by LTJG Greg Sabra, 

USCG. The steady force results were determined in a manner very similar to the method 
used to compute steady force test results for MIT tests. The method for computing MIT 
test results will be discussed in a later section. A complete discussion of LTJG Sabra’s 
code is contained in his thesis entitled ““Wave Effects on Underwater Vehicles in 
Shallow Water.” 

All of the MIT test conditions were listed in a common Excel file called “MIT Test 
Plan.xls”. This file contains separate worksheets for each of the many series of tests that 
were performed. These worksheets look very similar to the table contained in Appendix 
B, with the addition of two columns at the left to record date and time information for 
each experiment. The worksheets in the MIT Test Plan file were used by the MATLAB 
routines to determine which data files to analyze and what some of the test conditions 
were for each experiment. The test conditions obtained from the test plan were the water 
depth and submergence of the vehicle during the test. All other test conditions were 
extracted directly from the test data file. 

5.3.1. User Interface and Data File Management 

The user interface and the file management were performed by a MATLAB routine 
called “AutoanalyzeXls.m”. This file is contained as Appendix C. After the user starts 
this program, the user selects the series of experiments to be analyzed by entering the 
number corresponding to the desired series. All series that have been performed are 
listed, even those that are suspected to be of little value. After the test series has been 
selected, the program imports the list of experiments and the depth and submergence 
information. The program then calls other MATLAB routines to analyze the data files. 
For steady force tests, the analysis program is "AnalyzemodXlsSF.m". For all other tests, 
the analysis program is "AnalyzemodXls.m". 

The analyses were performed using the data files recorded by the notebook computer 
in Excel format. An example of a data file is included as Appendix D. 

5.3.2. Steady Force Data Analysis 

Steady force tests at MIT involved towing the vehicle down the tank with a steady 
angle of yaw or pitch. These tests were analyzed using "AnalyzemodXlsSF.m", 
contained in Appendix E. "AnalyzemodXlsSF.m" starts by importing the data file 
identified by “AutoanalyzeXls.m”. The program determines the ordered parameters and 
the date and time at which the experiment occurred. Then, the file eliminates the first 1.2 
seconds of data to allow for gantry acceleration and any data recorded after the gantry 
velocity returns to zero at the end of the test. The remaining position data is converted to 
the MKS system using a conversion factor. The remaining force and moment data 
undergo a more detailed analysis. 

The voltage output of the transducer is converted to forces and moments using 

F = - ^^ -- (14) 

GainxV^^^xSxlQ-^ 
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where 


F is the calculated force or moment 
Vout is the output voltage recorded by the computer 
Gain is the gain of that channel in the amplifier 
Vcxc is the excitation voltage of the channel, and 
S is the sensitivity of the channel. 

The calculation of equation (14) is performed using matrices so that the effects of cross¬ 
talk in the transducer can be accounted for. 

The mean force is calculated by taking the mean of the forces measured in the data 
interval and shifting the origin of the mean force from the origin of the load cell to the 
origin of the vessel coordinate system, defined to be at the midships of the vessel. The 
origin shift was done using 


^midships ^transducer ^^iransducer 

N =N +Yx 

midships transducer transducer 


(15) 


The origin-shifted mean forces and moments were written to a common output file that 
contained the mean force and moment data for all analyzed experiments. This output file 
and explanatory notes are included as Appendix F. 


5.3.3. Analysis of Experiments Involving Unsteady Motion 

The analysis of experiments involving unsteady motion was performed using 
“AnalyzemodXls.m” called by “AutoanalyzeXls.m”. The code is included as Appendix 
G. This is the most complicated of the codes used for this research and, as a result, is the 
most heavily commented. 

The first section of the code identifies and defines most of the variables used in the 
code. Next, the code initializes by reading the data file and gathering some basic 
information about the parameters of the experiment. The last four lines of the data file 
contain information about the ordered frequency and amplitudes of oscillation as well as 
sample frequency, velocity, and travel duration and distance. Then, the actual sample 
frequency is calculated by taking the inverse of the average interval between data points 
according to 

fsample 

mean(At) 

The ordered sample frequency was always 25 Hz, but for a certain period of time during 
the research errors in the control software resulted in data being taken at other 
frequencies. 

The code drops the first 1.2 seconds of data to allow for the acceleration of the gantry. 
The time interval to drop was chosen short enough to allow sufficient time remaining to 
have several periods of oscillation remaining but long enough to remove the majority of 
the acceleration period. The code also drops data recorded after the vessel completed its 
travel along the tank. This was necessary because the control software continued to 
collect data until the ordered time period of the experiment was completed, whether or 
not the travel distance had been accomplished. Setting the time period of travel too short 
resulted in sudden stops of the gantry causing large accelerations on both the vehicle and 
the gantry system. Ordered durations were made longer than absolutely necessary to 
prevent this mechanical shock to the system and prolong the life of the apparatus. 




The analysis code determines the frequency of oscillation by checking the input 
parameters in the data file to determine the ordered frequency of oscillation. This 
information is used to ensure that the data to be analyzed consisted of an integer number 
of wavelengths of the oscillation. This feature was absolutely necessary to get highly 
accurate results from the Fourier analysis that takes place later in the program. The 
period of a cycle is given by 

period =-^- 

frequency 

The duration of data recorded was found by taking the difference in time between the 
first and last remaining data points. The number of periods recorded is 

, - . , duration 

number of periods =-- 

period 

The number of data points retained for Fourier analysis is found by rounding down to the 
next integer the product of sample frequency, period, and the number of periods 
according to 

# of data points = round* period * y7oor(number of periods)) 

“round” is a MATLAB function that rounds the element to the nearest integer, “floor” is 
a MATLAB function that round the element to the next lower integer. 

Once an integer number of data points is established it the mean force and moments 
are subtracted from all measured forces and moments in order to remove steady effects. 

Next, the voltages from the transducer are converted to forces and moments using 
equation (14) and the numeric position data is converted to metric system position data 
using known relationships between the controller data and gantry motion. 

The force, moment and location data are conditioned by the program in preparation 
for the Fourier analysis. With the position data in metric format, the program translates 
the position data from the location of the strut to the vehicle midships. For linear motion, 
the motion of the strut forward of midships represents the motion of midships. For 
angular motion, this is not the case. The effect of angular motion on the x,y,and z 
position of midships is calculated by 

= distance sin 



where 


^midships ^strut ^strut Si^ 


^ 9 — 

180 j 


7 — — y 

^midships ^ strut strut 


1-COS 6 




TT 

m. 


(17) 


9 is the pitch angle 
If/ is the yaw angle 

distance is the distance along the x axis from the strut to midships. 

Lstrut is the length of the strut arm from its pivot point to the vehicle. 

The forces and moments are shifted from having their origin at the transducer to having 
their origin at midships using equation (15). Also, the data is interpolated into even 
intervals of exactly 0.4 seconds. The mean value of position for each channel except the 



X position is subtracted to remove any bias in the position data. Then, the data matrix is 
padded with zeros to obtain exactly 2048 data points. 

The most precise Fourier transformation requires the data to have an integer number 
of periods of the waveform and the frequency of signal to be analyzed must be a multiple 
of the fundamental frequency of the data sample. For an interval of 2048 data points 
being sampled at 25 Hz, the fundamental frequency is 

The frequencies of oscillation used for this research are 0.402831,0.79346, and 1.19629 
Hz representing 33, 65 and 98 times the fundamental frequency of the analysis. The 
sample frequency is assured by interpolating the data into exact time intervals between 
data points. As a result, the force, moment, and position data relating to oscillation are 
readily extracted using Fourier analysis. 

The first step in the Fourier analysis is to begin to build the data matrix by 
constructing the frequency column. The first row is assigned a frequency of zero Hz and 
each successive row is assigned a frequency of the row number multiplied by the 
fundamental frequency. The fast Fourier transformation is applied to the force, moment, 
and position data. To compensate for the zero padding added earlier, the value of each of 
the Fourier coefficients is multiplied by the ratio of the padded size to the unpadded size. 

After the Fourier transformation occurs, each coefficient has both a magnitude and 
phase associated with it. The phase of each of the coefficients is changed to make it 
relative to the phase of the motion that produced the force. This is done by multiplying 
every coefficient by where (p is the phase angle of the driving motion at that 
frequency. 

Low pass filters were installed on all of the force data collection channels to reduce 
the effects of electronic noise in the system. The effects of these filters is removed from 
each channel of force and moment using 

'nunfiUered ~ Vfiltered^ ^ ^ 

where the amplitude of the signal with filtering effects removed 

'Hfiltered ampHtudc of the signal after filtering 

a) is the frequency of oscillation 
C is the capacitance of the filter, and 
R is the resistance of the filter. 

The effects of filtering at the frequencies of oscillation altered the magnitude by less than 
one percent and the phase angle by less than five degrees. 

The program determines the frequency of motion closest to the ordered frequency of 
oscillation for all motions. The actual amplitudes of motion and forces are converted 
back into the time domain from the frequency domain using 


7(0 =. 


7(<y) 


where r]{t) is the amplitude in the time domain 

T]{co) is the amplitude in the frequency domain, and 
2048 is the number of data points used in the analysis. 
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The actual amplitudes of motion are used to calculate the inertial forces experienced 
by the vehicle using equation (2). The inertial forces are then subtracted from the 
measured forces and moments to leave only the hydrostatic forces and moments 
remaining. The mass, inertia, and center of gravity terms used to calculate the inertial 
forces and moments were derived by performing Oscillation tests in air. The model was 
filled with water to ensure that the effective mass of the vehicle was the same as it would 
be if the vehicle were in water. This guaranteed that the only significant forces measured 
were the inertial forces. The validity of the process was checked by performing 
oscillation tests in air and subtracting the calculated inertial components. The results of 
these inertial calculation checks are included as Appendix H. 

At this point in the program all of the hydrodynamic forces and moments at the 
frequencies of oscillation and their phases relative to the driving motion are determined. 
The results are conditioned in order to have positive amplitudes and have magnitudes of 
the phase angles less than 180 degrees. 

The program has also determined the forces and moments at twice and three times the 
oscillation frequency. In all cases the forces and moments at multiples of the oscillation 
frequency have magnitudes of approximately 10% or less of the forces and moments at 
the oscillation frequency. This indicated the response of the complete hydrodynamic 
system is linear and that non-linear forces and moments can safely be neglected in 
analyzing vehicle dynamics or in designing control systems. This also indicates that 
problems such as hysteresis in the load cell o-rings were unlikely. 

The results are written to an output file and contain all of the test information as well 
as the actual amplitudes and frequencies of oscillation and the amplitudes and frequencies 
of all six forces and moments. The phase angle between the force or motion and the 
driving is also listed. The output also includes the frequency, amplitude, and phase 
information for the second and third harmonics forces and moments. A partial example 
of the output file is contained in Appendix I. A complete output file consists of one row 
of 103 columns for each file analyzed. 

5.3.4. Curve Fitting the Experimental Results 

The following sections describe how the force and moment results were analyzed and 
present the resulting maneuvering coefficients. The effects of submergence and speed on 
the maneuvering coefficients will be analyzed. Part of that analysis will involve curve 
fitting the data to determine the functional relationships involved. The program used to 
do the curve fitting is a MATLAB routine generated by the author called 
“CoeffSolver.m”. T^e code is included as Appendix J. 

“CoeffSolver.m” uses user-coded values of the coefficients and the test parameters of 
Length/Submergence and Froude Number to perform a least squares regression of the 
data. When sufficient test data is present a second order equation is derived. When there 
is not sufficient data for a second order equation a first order equation is used. 

In general, the set of equations is of the form 

a, ,x, + 0,2^2 + 01,3X3 = b, 

^21^1 ^22^2 ~ ^2 

CljjXj + ^^32^2 ^33^3 ~ ^3 

which can be written in matrix form 






Ax = B 

where 

A is an m X n matrix 
X is an array of size m 
B is an array of size m. 

Matrix A represents the known parameters of the equations such as the speed and 

submergence of the vehicle for the experiment. Also, the terms aii usually equal 1 to 

allow for some constant to be built into the equation. The B array represents the 

measured or calculated data points. The array x represents the coefficients that are 

calculated to best represent the data using a least squares linear regression. The least 

squares regression calculates the elements of array x that will minimize the sum of the 

squares of the errors between the predicted and the calculated values. This manipulation 

is performed easily in MATLAB using x = A\B. 

To measure the quality of fit, the program finds the root mean square of the difference 

between the predicted and the measured coefficients. The code uses 

_ „ (Pr edicted - Actual) 

Difference =- 

Actual 

and 


Differences^, 


norm {Difference) 
sqrt{# of elements) 


where the norm of the Difference array is the largest singular value in the Difference 


array. 

Tbe output of the Coefficient Solver Program is included as 


5.4. Analysis of Experimental Results 


5.4.1. Rudder and Stern Planes Effects 

The effects of the rudder and stem planes were examined by performing tests at the 
United States Naval Academy with the control surfaces at no angle and with the control 
surfaces deflected to seven degrees. The results were normalized by subtracting the lift 
and moment at zero degrees firom the lift and moment measured with control surface 
deflection to remove any imbalance in loading. The lift coefficient per degree of fin 
deflection was calculated according to 

where Cls is the lift coefficient per degree. Lift is the measured lift force, U is the vehicle 
forward velocity, A^n is the effective area of the rudder, and 5 is the control surface 
angle. In order to provide more data points from the same experiments, the measured 
moments were converted into pseudo-lift forces by dividing the moments by the length of 
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the control surface moment arm, 0.7 m. These were then converted into coefficients 
using equation (21). 

For high aspect ratio wings, with aspect ratios greater than five, the theoretical 
approximation found by Hoemer’^ is 


da 


10 + 


20 


( 22 ) 


where AR^ is the aspect ratio found by 


^R. 


^ Spard ^ 

Area 


(23) 


The area was estimated by combining the calculated area of each of the fins with the 
estimated effective area provided by the body between the fins. Figure 9 shows the 
control surface geometry. Using an area of 0.02 m^ and span of 0.254 m, the aspect ratio 
is 3.23, and Cl is 0.0618/Degree. 



Figure 9. Control Surface Geometry 


For aspect ratios between three and five Hoemer recommends using 


'La 


da 


10 + 


10 
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(24) 


AR] .4R 

which yields Cl of 0.053/Degree.*'' 

Figure 10 shows the measured rudder lift coefficients plotted against the ratio of body 
length to submergence. The figure also shows the linear approximation to the data and 
the theoretical value for intermediate aspect ratio fins in an infinite fluid. The data for 
Figure 10 are presented in Table 3. 

The linear approximation to the data is 


= -0.0013 +0.0518 


Submergence 


(25) 


The value of the linear approximation at deep submergence is 0.0518/Degree, very near 
the theoretical value of 0.053/Degree calculated using equation (24). 
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Figure 10. Rudder Lift Coefficient for Various Values of Length/Submergence 


Table 3. Rudder Lift Coefficient Data 



Test 

Speed 

L/Subm. 

Rudder 

Angle 

Rudder 

Lift 

Normalized 
Rudder Lift 

Rudder 
Lift Coeff 

Yaw 

Moment 

Normalized 

Yaw 

Moment 

Yaw 

Moment 

Coeff 


m/s 


Degrees 

N 

N 

1/Degree 

N-m 

N-m 

1/Degree 

Measured Lift 

N1 

2.06 

2.07 

0 

-0.915 

0.000 


0.333 

0 


N6 

0.515 

4.14 

0 

-0.259 

0.000 

0.302 

0.000 

N7 

1.03 

4.14 

0 

-1.729 

0.000 

1.336 

0.000 

N8 

2.06 

4.14 

0 

-7.349 

0.000 

6.402 

0.000 

N9 

2.06 

4.14 

7 

12.267 

19.616 

0.0662 


-17.316 

0.0825 

N15 

0.515 

8.26 

0 

0.280 

0.000 


0.196 

0.000 


N16 

2.06 

8.26 

0 

2.810 

0.000 

1.512 

0.000 

N22 

2.06 

8.26 

7 

8.001 

5.191 

0.0175 

-10.594 

-12.106 

0.0577 


2.06 

8.26 

7 

7.922 

5.112 

0.0172 

-13.083 

-14.595 

0.0695 


N9m 

2.06 

4.14 

7 


-7.832 

0.0264 


Lift Calculated 
from Moment 

N22m 

2.06 

8.26 

7 

-17.294 

0.0583 

N24m 

2.06 

8.26 

7 

-20.850 

0.0703 


The effects of stem planes were determined in a manner similar to those of the rudder 
with similar results. Figure 11 shows the measured Stem Planes Lift Coefficients plotted 
against the ratio of body length to submergence. Figure 11 also shows the linear 
approximation to the data and the calculated theoretical value. The data for Figure 11 are 
included as Table 4. 
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Figure 11. Stern Planes Lift Coefficient for Various Values of Length/Submergence 


Table 4. Stern Planes Lift Coefficient Data 



Test 

Speed 

L/Subm. 

Stem 

Planes 

Angle 

Stern 

Planes 

Lift 

Normalized 

Stem 

Planes Lift 

stern 
Planes 
Lift Coeff 

Pitch 

Moment 

Normalized 

Pitch 

Moment 

Pitch 

Moment 

Coeff 


m/s 




N 

1/Degree 

N-m 

N-m 

1/Degree 

Measured Lift 

am 



0 

2.850 

0.000 


-3.227 


wm 

EE3 


7 

19.131 

16.281 

0.0549 

4.982 


0.0391 



4.14 

0 

EB3 

0.000 


-0.157 


am 


4.14 

0 

mm 


-0.778 



■lyi 

0 



-6.761 

N10 



-7 



0.0601 

-8.460 

-1.699 

0.0081 



B;IJJ 

0 

lilSBi 

BMi!iM 


-0.162 


jsm 

S3 


0 

SB3 


-8.868 

MB 

S3 


-7 

BIWtU! 

aEBEB 

0.0510 


Mlt;b!iB 

lifiSCT 

MB 

S3 


-7 


MIEB 

0.0712 



BisnsEi 

Lift Calculated 
from Moment 

Bl 

2.06 

2.07 

7 

1 


0.0396 


N20m 

2.06 

8.26 

-7 


0.0563 

N24m 

2.06 

8.26 

-7 


0.0712 


These results seem to show that the rudder has a slightly reduced effect when near the 
surface and the stem planes have a slightly greater effect when near the surface. In both 
cases, the linear approximation of the data closely approaches the theoretical value for a 
submerged body in an infinite fluid. The opposite slopes of the linear approximations are 
not explainable at this time. The small variations shown in the coefficients over the range 
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of submergences monitored show that the effect of submergence is very small for 
submergences greater than 10 % of body length. 


5.4.2. Sway Force and Yaw Moment Due to Sway Motion 

Experiments to determine the effects of sway motion were performed in various 
combinations of submergence, velocity, frequency and amplitude of oscillation. Table 5 
contains the test conditions and key results. 


Table 5. Test Conditions and Results for Sway Force and Yaw Moment Due to Sway Motion 


1 Sway Motion 

Hydrodynamic Sway Force, Y | 



Frequency 


Amplitude 


■asm 

Yvdot 

Yvdot’ 

■nasn 

Yv 

Yv’ 



Hz 

m 

N 


N 

_kg_ 


N 

msism 


0.543 


0.40283 

0.1 

2.217 




-0.0175 

-1.1944 

-4.7044 


0.252 

■lE&S 

0.40283 

0.1 

2.641 





-1.6937 



0.398 


0.79346 

0.1 

8.644 


mxmm 

-2.7725 


-5.0439 

«]tSl6WCM 


0.543 

■JliisUi 

0.40283 

0.1 

2.677 

-38.1 


-3.2812 


-1.6517 



0.252 


0.40283 

0.1 

■EQ9H 



-2.9156 

-0.0175 

-1.1841 

-4.6652 



1 Sway Motion 

Hydrodynamic Yaw Moment, N 1 



Frequency 

Qsnioss 


Phase 

WESsim 

Nvdot 

Nvdot’ 


■C9H 

Nv’ 

m 

m/s 

Hz 

m 

N-m 


■GSH 

_kg_ 


■dsn 

■SBI 


0.543 


0.40283 

0.1 

0.085 

-45.4 




EEm 



0.252 

IKiMiM 

0.40283 

0.1 

0.351 


0.0513 



-0.34714 

-1.36874 


0.398 


0.79346 

0.1 

0.352 

-55.8 



giTiimna 



BiPliEtni 

0.543 

Hblslsl 

0.40283 

0.1 

0.371 

-80.1 





-1.44166 


0.252 

lilcgRM 

0.40283 

0.1 

0.085 







B.M.UEI 


The hydrodynamic forces and moments are determined by subtracting the inertial 
forces and moments from the measured forces and moments for each experiment. For 
pure sway motion, the applicable equations of motion from equation (9) are 

Y = Yv + Y,y 


(26) 


(27) 


N = N^v+N,v 

The notation of complex equations is used to separate the components of the measured 
force and moment into the component related to velocity, a damping component, and the 
component related to acceleration, an added mass component. This is done using 

T = Re(T) + /• Im(y) = Y cos((f>) + iY sin(^z5) 

N = Re( Y) + i lm(N) = N cos(^) + iN sin(^) 

and understanding that 

Re(r) = rcos(^^) = y,v 
Im(y) = ysin(<z)) = y„v 
Re(7V) = Arcos(^^) = y,v 
Im(Y) = 7Vsin(^i) = y„v 


(28) 


^ where (j) is the phase angle taken by subtracting the phase angle of the force from the 

phase angle of the driving motion. 

Sinusoidal motions can be described by 

r] = rie‘“‘ (29) 

where 7 is a time varying position 

# 7 is the amplitude of the sinusoidal oscillation, and 

(o is the frequency of oscillation. 
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The velocity component is 


(30) 


^^coAe“ 

dt 


and the acceleration component is 


d^rj 

IF 


= -Q)^Ae“ 


(31) 


In the standard method of notation the exponential term is understood and velocity is 
represented by coA and the acceleration term is represented by ~co^A. 

Therefore, 

V = CO A 


(32) 


V = -co^A 

By combining equations (28) and (32), we obtain the necessary equations to compute the 
hydrodynamic forces and moments due to sway according to , 

Im(y) 


N,= 


(oA 

Re(y) 

-co^A 

ImjN) 

coA 

Re(A^) 

-co^A 


(33) 


The results are non-dimensionalized according to the method of Section 4. 

5.4.3. Heave Force and Pitch Moment Due to Heave Motion 

Experiments to determine the effects of heave motion were also performed in various 
combinations of submergence, velocity, frequency and amplitude of oscillation. Table 6 
contains the test conditions and key results. 

Table 6. Test Conditions and Results for Heave Force and Pitch Moment Due to Heave Motion 


t Heave Motion 

Hydrodynamic Heave Force, Z I 


OSESG9 


Amplitude 

Amplitude 




Zwdot' 

}m(Z) 

Zw 

Zw' 

m 


Hz 

m 

N 


N 

_kg 


N 



0.488 


0.40283 

0.1 

2.406 


2.0452 


-0.0192 

-1.2681 


BiliSga 

WKEEBKKk 


0.40283 

0.1 


mKorm 

WBSSEMt 

-3.8243 

-0.0230 


-6.9812 


0.398 


0.79346 

0.1 



6.9747 


-0.0167 


BEssa 


0.488 

IgiliM 

0.40283 

0.1 

2.687 



-3.4543 


-1.5100 

-5.9418 

emei 

0.272 

■filikkli 

0.40283 

0.1 

2-633 

-28.9 

2.3054 



-1.2727 

■KrtYriiai 

gili^ 


1 Heave Motion 

Hydrodynamic Pitch Moment, M | 

Submergence 






■i^mi 

Mwdot 

Mwdot’ 



Mw* 

m 

nfi/s 

Hz 

m 



N-m 

i^g 


■SBuH 

■UiH 


niEEnisii 


0.40283 

EEnH 


■FSIililil 

■BSIEEISi 

0.32703 

mmmm 

E^iEI 

1.02217 


0.272 


0.40283 

0.1 

0.395 

■mBiaa 

-0.1472 

0.23097 



1.45422 

iimi];r£i 

0.398 


0.79346 

0.1 

0.966 


-0.8223 

0.32843 

MiHthWA-M 

■tuatiasiB 

1.00733 


0.488 

MIiW 

0.40283 

0.1 



■E&SISI 




1.75554 


0.272 

0.333 

0.40283 

0.1 

0.293 


-0.2372 

0.37125 

0.00322 

0.17110 

■ilHririrfil 

■miKWci 


The analysis of the results in heave motion follows t 


he same train of reasoning as 


described previously for sway motion. The equations for heave motion become 


Z = Z^w + Z> 

M = M^yv + M^w 


(34) 
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and 



Im(Z) 

coA 


Z... 


M,. 


M, ,= 


Re(Z) 

-(O^A 

Im(M) 

6)A 

Re(M) 


( 35 ) 


5.4.4. Discussion of Heave and Sway Motion Results 

The results of these experiments are in agreement with the expected results. Because 
the sway and heave motions with forward velocity create effective angles of attack of the 
body, the bow and stem both experience lift force opposed to v and w, therefore and 

Z^ are always negative. The terms and Z^ are always negative and have a magnitude 
approximately equal to the displacement of the vessel. and are usually negative. 


but can become positive if the rudder or stem planes are very large. and are 

usually relatively small quantities of uncertain sign.’* The results shown in Table 5 and 
Table 6 match the expected values. 

The results contained in Table 5 and Table 6 are graphically presented in Figure 12 
through Figure 27. The results have been made non-dimensional according to the method 
described in section 4. 
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Figure 12. YJ and Z^ vs L/Subm at Froude Number = 0.128 





Damping vs L/Subm at Froude Number = 0.128 
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Figure 13. YJ and ZjvsL/Subm at Froude Number = 0.128 


Added'Mass vs L/Subm at Froude Number = 0.383 



Figure 14. Y-' and Zj vs L/Subm at Froude Number = 0383 
















Damping vs L/Subm at Froude Number = 0.383 
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Figure 15. YJ and Zj vs L/Subm at Froude Number = 0383 


Added Mass vs Froude Number at L/Subm = 1.277 for 
Sway and 1.42 for Heave 
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Figure 16. and Z/ vs Froude Number at L/Subm = 1.277 for Sway and 1.42 for Heave 















Damping vs Froude Number at L/Subm = 1.277 for Sway 

and 1.42 for Heave 



Figure 17. YJ and Zj vs Froude Number at L/Subm = 1.277 for Sway and 1.42 for Heave 


Added Mass vs Froude Number at L/Subm = 2.749 for 
Sway and 2.547 for Heave 
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Figure 18. F- and Z^, 


VS Froude Number for L/Subm = 2.749 for Sway and 2.547 for Heave 

















Damping vs Froude Number at L/Subm = 2.749 for Sway 
and 2.547 for Heave 



Figure 19. Y' and Z^, vs Froude Number for L/Subm = 2.749 for Sway and 2.547 for Heave 


Figure 12 through Figure 15 show that the effect of submergence on the direct added 
mass and damping terms for sway and heave motion is nearly negligible. Figure 16 
through Figure 19 show that there is a significant effect of speed on these terms. These 
figures also show that the direct terms in sway and heave behave very similarly. The 
cross-term added mass and damping coefficients shown in Figure 20 through Figure 27 
have different patterns of behavior 










Figure 21. NJ and MJ vs L/Subm for Froude Number = 0.128 
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Figure 22. and vs L/Subm for Froude Number = 0.383 
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Figure 23. N' and M' vs L/Subm for Froude Number = 0.383 
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Cross Term Added Mass vs Froude Number for L/Subm = 
2.749 for Sway and 2.574 for Heave 
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Figure 26. N’ and vs Froude Number for L/Subm = 2.749 for Sway and 2.574 for Heave 



The cross-terms due to sway and heave motion are expected to be very similar to each 
other and have similar dependencies, however Figure 20 through Figure 27 show that this 
may not be the case. Examination of the measured moments listed in Table 5 and Table 6 
show that the measured yaw moments due to sway force range from 0.085 to 0.352 N-m 
and the measured pitch moments due to heave force range from 0.293 to 0.966 N-m. The 













transducer is rated for up to 28 N-m in yaw and pitch. This means that the maximum 
measured moments represent less than 4% of the range of the transducer. The minimum 
moment represents only 0.3% of the range of the transducer. The abnormal behavior of 
the moment coefficients may be due to insufficient moment being present to properly 
deflect the transducer. 

The effects of submergence and speed were determined numerically by performing a 
least squares fit to the data using linear regression. This procedure is described in section 
5.3.4. The resulting equations for the coefficients are: 

Y' = -0.081458 + 0.001774—— + 0.121754Fr 

Submergence 

Y' = -0.016345 + 0.000061—-o.007224Fr 

'' Submergence 

Z; =-0.086685 + 0.00091—-+ 0.140751Fr 

Submergence 

z/ = -0.013633 - 0.002679—-0.005671Fr 

Submergence 

M„; = 0.023224 - 0.002943—-0.022359Fr 

Submergence 

M / = 0.002825 + 0.000594 —-0.006407Fr 

Submergence 

N' = -0.00207 + 0.000093—-0.016489Fr 

'' Submergence 

N’ = -0.00089 + 0.000037 — -EEJ I -+ 0.000172Fr 

Submergence (36J 

The quality of the fit for these equations is expressed in terms of the root-mean-square 

value of the percent difference between the predicted and the empirical coefficients. 

Table 7 contains the values for the quality of fit. 


Table 7. Quality of Fit for Heave and Sway Motion Coefficients 


Coefficient 

Fit 


15% 


4% 


17% 


9% 

< 

19% 

m: 

27% 

K 

9% 

n; 

9% 













5.4.5. Forces and Moments Due to Body Angle 

The forces and moments due to body angle were determined using data from 
experiments performed on the small scale model at MIT. Appendix L contains the list of 
experiments performed. 

The equations for the coefficients due to body angle are'^ 

Y 


Y..=- 


-^pU^D^a 


Z .=■ 


M... =- 


I pU^D^a 
M 

—pU^D^La 


(37) 


where 


N... 


N 


— pU^D^La 


D is the diameter of the vehicle 
L is the length of the vehicle 
a is the angle of attack. 


The results of these experiments are included as Appendix M. The effect of the angles 
is assumed to be linear for small angles, so the coefficients were calculated for each 
experiment and then averaged for each combination of speed and submergence. Pitch 
and Yaw angle coefficients were averaged separately. Any data point with a calculated 
coefficient more than 1.15 standard deviations from the mean was removed from 
consideration as unreliable data. 1.15 standard deviations was chosen as the 
discrimination point in order to remove the clearly bad data while retaining as much of 
the possibly good data as possible. The mean coefficients for each combination of 
submergence and velocity are shown in Table 8. Figure 28 through Figure 31 illustrate 
the dependency of the restoring forces on submergence and speed. 

It is very important to note at this point that the magnitudes of the yaw and pitch 
moments measured during these tests are very small, on the order of less than 1% of the 
capacity of the load cell. The data is analyzed and presented here, but further work is 
required to determine the accuracy of the equations with a more appropriate transducer. 
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Table 8. Mean Coefficients for Force and Moment Due to Body Angle at Various Submergences and 
Velocities 


Submergence 

Velocity 

Yuv 

m 

m/s 


0.252 

0.75 

-5.36887 

0.398 

0.75 

-4.85345 

0.543 

0.75 

-2.32929 

0.252 

1 

-3.51475 

0.398 

1 

-2.8487 

0.543 

1 

-1.61492 


Zuw 

Submergence 

Velocity 

Zuw 

m 

m/s 


0.252 

0.75 

1.023955 

0.398 

0.75 

2.269136 

0.543 

0.75 

1.141877 

0.252 

1 

1.220249 

0.398 

1 

2.612357 

0.543 

1 

1.801033 


Submergence 

Velocity 

m 

m/s 

0.252 

0.75 

0.398 

0.75 

0.543 

0.75 

0.252 

1 

0.398 

1 

0.543 

1 


-1.02061 


-0.99671 


-0.68264 


-0.80028 


-0.73773 


-0.58644 


Submergence 

Velocity 

m 

m/s 

0.252 

0.75 

0.398 

0.75 

0.543 

0.75 

0.252 

1 

0.398 

1 

0.543 

1 


-0.626 


-0.59447 


-0.71149 


-0.4012 


-0.61916 


-0.59916 
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USubm 



'Fr = 0-290 - -■ — ’Fp=0.386 


Figure 29. Zuw as a Function of Submergence for Two Speeds 



L/Subm 


^ -*Fr= 0.290 — ’Fr=0.386 


Figure 30. Nuv as a Function of Submergence for Two Speeds 










Figure 31. Muw as a Function of Submergence for Two Speeds 

The data from the restoring force and moment experiments was put into a linear 
regression model to determine the dependency of those forces and moments on 
submergence and speed. The resulting equations are 

r, =3.035943-11.232476—— + 2.399695 f—-1 +15.79518 

Submergence Submergence j 


+ 15.795188Fr 


= -7.723764 + 9.156976 -2.365368 

Submergence 

=0.020017-1.452965—— + 0.317964 

Submergence 


Length 

^ Submergence 
Length 

Submergence j 


+ 4.182974Fr 


+ 1.987686Fr 


= -1.168727 + 0.127681 —-0.007568 

Submergence 


Length 

Submergence 


+ 1.079278Fr 
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5.4.6. Sway Force and Yaw Moment Due to Yaw Motion 

Experiments were performed with several combinations of submergence, velocity, and 
frequency of oscillation to determine the effect of yaw motion on sway force and yaw 
moment. The experiments performed are listed in Appendix B. The method of analysis 
was similar to that used for the sway and heave motion tests, but in this case there is 
another term due to the angle of the body as it moves forward and oscillates in yaw. 

The equations are 


and 


N = N,r + N,r + N^{-\pU^D''La) 
Im(F) 


(39) 


F = 


F.= 


AT = 




coA 

Re(F)-F„,(-iy9l7^DV) 

-gP'A 

Im(Y) 

(oA 

^t{N)-Y^{-\pU^D^La) 

-coi^A 


(40) 


The restoring force is a real force and must be subtracted from the measured real 
hydrodynamic force to calculate the force due to the yaw motion. 

The test conditions and results are presented in Table 10. 

Table 10 Test Conditions and Results for Sway Force and Yaw Moment Due to Yaw Motion 


t Yaw Motion 

Hydrodynamic Sway Force, Y I 

Submergence 






■UMOBi 

umsm 

Yrdot 

Yrdot’ 

Im(Y) 

Yr 

Yr* 

m 


Hz 


N 


N 

N 

kg 


N 

kg/s 


0.252 



10 

»l;gf 

ggngnrii 

-1.0402 


0.1543 

0.0013 

1.5192 

1.1650 

0.0210 

0.252 



10 


mitrmm 

-0.3928 

-2.4287 

0.2462 

0.0021 

3.1095 

2.3690 

0.0142 

0.543 



10 


mesem 


MignaM 


WKimsm 

MtgKIB 



0.543 



10 

1.760 

Igiailil 

-0.9994 

-1.3294 

0.1359 



IKKliciil 


0.398 



10 

1.761 

KlUiiia 

0,9565 

-0.3364 

0.3010 

■ililiha 




0.543 



10 


HiUCMil 


-1.0481 

0.2425 

0.0021 


Ulicfel 

■iliHiZl 


Yaw Motion 

Hydrodynamic Yaw Moment, N | 

Submergence 

Velocity 

Frequency 

Amplitude 

Amplitude 

Phase 

Re{N) 

Re(N>-Nuv 

Nrdot 



■m 

■cm 

m 

IBu!SiH 

Hz 


HE3C31I1 


N-m 

N 

kg 


N-m 

■ss 


0.252 


■QQ9I 

10 





-0.0650 


ggRRgg 


faikiiiKkii 


■KiWiM 

1.196 

10 

1.114 

-34.2 






gtierir^rti 


0.543 



10 

■iKW 

-21.2 

WKSESBMM 




giliCTZiH 



0.543 

■tlcjckil 

1.196 

10 

■ig;W 


0.7008 


-0.0670 


ESSEEI 


BIIIIIKUM 

0.398 



10 


-22.7 

0.3229 




BilkkfiVJ 



0.543 



10 


-27.8 

■ESEEQHi 




Ei&EILSI 




5.4.7. Heave Force and Pitch Moment Due to Pitch Motion 

The analysis for heave force and pitch moment due to pitch motion closely mirrors the 
analysis for sway force and yaw moment due to yaw motion. For pitch motion, the 
equations are 


Z = Z^q + Z.q + Z^S-\pU^D^a) 
M=Ma + MA + M^S-\pU^-D^La) 


(41) 


and 
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Im(Z) 

^ coA 

" -a)^A 

M 

coA 


(42) 


Rc(M)-M^J-^pU^D^La) 

’ -6>^y4 

The test conditions and results are presented in Table 11. During these tests special 
consideration was given to removing the effects of surge motion caused by the long pitch 
arm of the test apparatus. In order to remove these effects, the test apparatus was 
oscillated in surge at the same time pitch oscillations occurred. The surge oscillations 
were 180° out of phase with the pitch oscillations and of such a magnitude as to cancel 
the surge due to pitch. 


Table 11. Test Conditions and Results for Heave Force and Pitch Moment due to Pitch Motion 


1 Pitch Motion 

Hydrodynamic Heave Force, Z I 

Submergence 






Re(Z) 




mSimm 

zq 

ZqV 

m 


■BIOH 


WKEMM 


N 

N 



HBSHi 

■!Mm 


0.252 


1.19629 

9 

■clMtM 






mmm¥M 



0.252 


WEEmm 

9 




-0.4326 


■am 




0.543 

■•mi 


9 


KlXlilii 


WKES3Mtl 

giKtdcgl 

gtititiga 



BtygkRg 

0.543 

■iIckM 

WEMESm 

9 






MwhMM 



giltkglil 

0.252 

»mnif 


9 


MSfiTiCT 

0.1579 

■WiyLhi.M 

mn^mm 




BtliliEa 


I Pitch Motion 

Hydrodynamic Pitch Moment, M I 

Submergence 

■Miifggiga 

OSffi55!l 

ESSI^SISSSI 








■CSuB 

mSSm 

m 


Hz 

ai!5«ff53M 

■dSB 


N-m 

H 

_Isa_: 


N-m 

UB 


0.252 


■SEBI 

9 



0.5470 

0.5749 


jgilifiliUM 

Eilc?iRgl 

BilggZtn 

giitTitnai 

0.252 



9 

KEBI 


IHQQSSi 

WKSESSMi 

gtyjcliKa 


EiHWiliia 

giimiTiii 


0.543 

■flgklcM 

1.196 

9 


^msi 

■CESiB 




BtgirtHil 

giRIAIigl 


0.543 

■iIckkM 


9 




heeebi 

■riiffiTga 


gqggga 

giMEftll 


0.252 



9 


g|JiIi!il 

0.1062 

MiKWikl 


giTiTiMJ 

mhmbi 

giIiW;k>i 

gililiWi!:! 


5.4.8. Discussion of Yaw and Pitch Motion Results 

The results of yaw and pitch motion testing are graphically illustrated in Figure 32 
through Figure 47. Both submergence and speed were found to be significant factors. 
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Added Mass vs L/Subm at Froude Number = 0.128 
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Figure 32. and M^ vs L/Subm at Froude Number = 0.128 


Damping vs L/Subm. at Froude Number = 0.128 


0.00000 

-0.0010Q. 

- 0.00200 

-0.00300 

-0.00400 

-0.00500 

-0.00600 

-0.00700 

-0.00800 

-0.00900 

- 0.01000 


000 - 


-4:500- 


-^. 000 - 


2.500 


-O.0018X- 0.0035 


y = -0.0007x-~'O; 


3:000 



+ Nr' 

X Mq' 

-Linear (Nr') 

-- Linear (Mg') 


Length/Submergence 


Figure 33. and M’ vs L/Subm at Froude Number = 0.128 










Added Mass vs L/Subm at Froude Number = 0.383 
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Figure 34. N-' and M^ vs L/Subm at Froude Number = 0383 


Damping vs L/Subm. at Froude Number = 0.383 
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Figure 35. N' and vs L/Subm at Froude Number = 0.383 













Added Mass vs Froude Number for L/Subm = 1.277 
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Figure 36. and vs Froude Number at L/Subm= 1.277 
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Figure 37. and M vs Froude Number at L/Subm= 1.277 













Added Mass vs Froude Number for L/Subm = 2.749 
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Figure 38. N-' and vs Froude Number at L/Subm= 2.749 


Damping vs Froude Number for L/Subm = 2.749 
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Figure 39. and vs Froude Number at L/Subm= 2.749 


















Figure41. Y* and vs L/Subm at Froude Number ==0.128 
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Figure 43. YJ and Z^' vs L/Subm at Froude Number = 0.383 






















Cross Term Added Mass vs Froude Number at L/Subm = 

2.749 


0.0050 


0.0040 + 


0.0030 


0.0020 

0.0010 


y=0.0202x-0.0041^ 
-—-^ 


0.0000 

- 0 . 0015 - 


= 0.0099x + 7E-05 


+ Yrdof 
X Zqdof 

Linear (Yrdof) 
— — Linear (Zqdof) 


.000 0.100 Wf.200 0.300 0.400 


- 0.0020 


Froude Number 


Figure 46. Y^. and Z- vs Froude Number at L/Subni= 2.749 
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Figure 47. YJ and Z ' vs Froude Number at L/Subm= 2.749 


The data from pitch and yaw motion was also analyzed by linear regression to yield 
equations for the coefficients. Those equations are 













V ' = 0.020906 + 0.000403 -0.017979Fr 

Submergence 

Y\= 0.002324 - 0.000087 -o.000632Fr 

Submergence 

Z ’ = -0.070199 + 0.003498—— + 0.133729Fr 
^ Submergence 

Z.' = -0.002666 - 0.000514—— + 0.020222Fr 
’ Submergence 

M ' = -0.010515 - 0.000701 — + 0.030777Fr 

^ Submergence 

M, ' = 0.00014 + 0.000106—-0.009854Fr 

^ Submergence 

' = -0.006013 - 0.001173—— + 0.012297Fr 

Submergence 

N, ' = -0.000786 - 0.000019 —^ _ 0.000582Fr 

Submergence 

Table 12 contains the values for the quality of fit. 

Table 12. Quality of Fit for Pitch and Yaw Motion Coefficients 


Coefficient 

Quality of Fit 

i;' 

7% 

n' 

46% 

z; 

0% 


0% 


0% 


0% 

< 

23% 

n; 

8% 


(43) 


6. Computational Analysis 

Numerical method results to accompany experimental results are very important. To 
that end the validated free surface linear code for surface ships, SWAN, was revised by 
one of its developers to include submerged objects. Unfortunately, the results from the 
submerged vehicle version of SWAN were too unreasonable to be either valid or 
included here. For example, in some conditions, the added mass was predicted to be 16 
times the actual displaced mass. That cannot be correct. There was not enough time to 
properly revise, test, and validate the new version of this numerical method. Zero 
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# 


♦ 





forward speed evaluations of added mass and damping were performed with the validated 
code WAMIT and those results were reasonable. However, the interest for this research 
is for an underwater vehicle with forward speed. 

7. Conclusion 

The analyses performed for this research have drawn upon a limited number of 
experiments with a great deal of uncertainty in the quality of data. The limited number of 
experiments is a result of very significant delays experienced in designing and testing the 
test apparatus. Although limited data exists that is known to be valid, a very large 
amount of data exists from earlier experiments. These earlier experiments may have had 
interaction between the transducer and the vehicle shell. Although the geometry of the 
model indicates no interference, there may have been interference due to shell distortion 
and corrosion products. Also, some of the uncertainty in the quality of early data results 
from the measurement of very small forces and moments that represent a very small 
fraction of the range of the transducer. The method of research and analysis is sound and, 
at the very least, provides a roadmap for conducting the research in the future. 

The results associated with the sway force due sway motion and heave force due to 
heave motion are very similar and indicate that the quality of that data is very good. The 
effects of variation in submergence and speed on those forces can also be expected to be 
good results. 

The proximity to the free surface has relatively little effect on forces due to motions or 
control surface deflections. This indicates that with minor control system alterations, 
operation in shallow water is probably feasible. Of course further work on proximity to 
the bottom is necessary. 

The causes for the variation in the results with Froude number and submergence are 
not well understood. There are various theories that discuss the interaction of the flow 
stream with the free surface and the bottom that may be useful, but true understanding of 
the results of these analyses must be saved for future work. 

8. Future Work 

In this thesis we have found the effects of shallow water and submergence on many of 
the dominant terms in the maneuvering equations. The reason for these variations is also 
yet to be determined. In the very near future the data from the experiments conducted 
early in this research and not known to be reliable will be analyzed. If that data is found 
to be reliable, the results reported here will be expanded and modified to include that data 
and provide a better description of the effects of submergence and speed. 

This work indicates that the proximity to the free surface has relatively little effect on 
forces due to motions or control surface deflections. This indicates that with minor 
control system alterations, operation in shallow water is probably feasible. Of course 
further work on proximity to the bottom is necessary. Much of this is available from the 
work of William Ramsey.'^ 

Additional future work will involve finding a numerical method that can accurately 
predict the phenomena studied here. The results of this and similar research will be used 
to validate that method for widespread use in predicting this kind of behavior. 
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Appendix A. Full Scale Experiments Performed at the 
United States Naval Academy to Determine the Effects of Body 


Angle and Control Surface Deflection 

NAVAL ACADEMY REMUS NO WAVE STEADY FORCE TEST RUNS 
_(depth is measured to top of vehicle at strut center) 


# 

• 

• 

# 

• 

Run 

# 

Rudder 

Angle 

(deg) 

Stern 

Plane 

Angle 

(deg) 

Pitch 

Angle 

(deg) 

Yaw 

Angle 

(deg) 

Speed 

(ft/s) 

Depth 

(meters) 

V/ 

sqrt(gL) 

depth/L 

N1 

0 

0 

0 

0 

6.8 

0.762 

0.526 

0.481 

N2 

0 

0 

0 

8 

3.4 

0.762 

0.263 

0.481 

N3 

0 

0 

8 

0 

6.8 

0.762 

0.526 

0.481 

N4 

0 

7 

0 

0 

6.8 

0.762 

0.526 

0.481 

N5 

7 

0 

0 

0 

3.4 

0.762 

0.263 

0.481 

N6 

0 

0 

0 

0 

1.69 

0.381 

0.131 

0.240 

N7 

0 

0 

0 

0 


0.381 

0.263 

0.240 

N8 

0 

0 

0 

0 

6.8 

0.381 

0.526 

0.240 

N9 

7 

0 

0 

0 

6.8 

0.381 

0.526 

0.240 

N10 

0 

7 

0 

0 

6.8 

0.381 

0.526 

0.240 

N11 

0 

0 

4 

0 

6.8 

0.381 

0.526 

0.240 

N12 

0 

0 

8 

0 

6.8 

0.381 

0.526 

0.240 

N13 

0 

0 

0 

4 

6.8 

0.381 

0.526 

0.240 

N14 

0 

0 

0 

8 

6.8 

0.381 

0.526 

0.240 

N15 

0 

0 

0 

0 

1.69 

0.1905 

0.131 

0.120 

N16 

0 

0 

0 

0 

6.8 

0.1905 

0.526 

0.120 

N17 

0 

0 

0 

8 

6.8 

0.1905 

0.526 

0.120 

N18 

0 

0 

8 

0 

6.8 

0.1905 

0.526 

0.120 

N19 

0 

0 

-8 

0 

6.8 

0.1905 

0.526 

0.120 

N20 

0 

7 

0 

0 

6.8 

0.1905 

0.526 

0.120 

N21 

0 

7 

-8 

0 

6.8 

0.1905 

0.526 

0.120 

N22 

7 

0 

0 

0 

6.8 

0.1905 

0.526 

0.120 

N23 

7 

0 

-8 

0 

6.8 

0.1905 

0.526 

0.120 

N24 

7 

7 

0 

0 

6.8 

0.1905 

0.526 

0.120 


62 
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Appendix B. Model Scale Experiments Performed at MIT to 
Determine the Forces and Moments Due to Unsteady Motion 





Velocity 


lOSClOOl 0.795 I 0.252 | _ | j j _|_ 1 _ j 10 j 0.40281 0.333 | 





































































































Appendix C. AutoanalyzeXls.m 

% AutoanalyzeXls.m 
% Erik Oiler, 2003 

% This program uses other programs to automatically, analyze all the data 
% files. 

% Operational Overview 

% 1. User starts this program to perform analysis of experimental data 
% files. 

% 2. User selects the series of experiments to be analyzed. 

% 3. The program opens an excel file called "MIT Test Plan.xls" and 
% imports the data from the worksheet for the selected series. The 
% matrix seriesnum contain the numerical data from the worksheet 
% and the matrix seriestext contains the text data from the 
% worksheet. 

% 4. The program runs the appropriate analysis program for each test 
% series. For steady force tests, the analysis program is 
% "AnalyzemodXlsSF.m". For all other tests, the analysis program 
% is "AnalyzemodXls.m". 

% 5. The called analysis program analyzes the raw data and writes the 
% will display that the data file does not exist. 

% 6. When all data files listed in the test plan have been analyzed, 

% file. 

% 7. The following files must be in the same directory: 

% "AutoanalyzeXls.m" 

% "AnalyzemodXlsSF.m" 

% "AnalyzemodXls.m" 

% "MIT Test Plan.xls" 

% All Data Files to be analyzed. 

% Initialize the workspace by clearing all variables and closing all 
% windows, 
close all; 
clear all; 


% Determine which set of tests to analyze. 

fprintfC 1: Horizontal Plane \n') 

fprintfC 2: Vertical Plane \n') 

f^rintfC 3: Pure Sway \n') 

fprintf(' 4: Pure Heave \n') 

fprintf(' 5: Pure Pitch \n') 

fj^rintfC 6: Pure Yaw \n') 

fprintfC 7: Mass Matrix \n') 




fprintfC 8: Steady Force \n') 
fprintfC 9: Inertial Calculation Checks \n') 
fprintfC 10: Miscellaneous \n') 
fprintfCl 1: Oscillation Tests \n') 

series = input('Which Test Series? \n'); 

% Select the worksheet in "MIT Test Plan.xls" and the output file based 
% upon the test series, 
switch series 
case 1 

worksheet = Hor Plane'; 
outputfile = 'HorPlaneOut.txt'; 
case 2 

worksheet = 'Compensating Vert Plane'; 
outputfile = 'VertPlaneOut.txt'; 
case 3 

worksheet = 'Pine Sway'; 
outputfile = 'SwayOut.txt'; 
case 4 

worksheet = 'Pure Heave'; 
outputfile = 'HeaveOut.txt'; 
case 5 

worksheet = 'Pure Pitch'; 
outputfile = 'PitchOut.txf; 
case 6 

worksheet ='Pure Yaw'; 
outputfile = 'Y awOut.txt'; 
case 7 

worksheet = 'Mass Matrix Tests'; 
outputfile ='MassTestsOuttxt'; 
case 8 

worksheet = 'Steady Force'; 
outputfile = 'SteadyForceOut.txt'; 
case 9 

worksheet = 'Mass Matrix Tests'; 
outputfile = 'InertiaCalcTestsOut.txt'; 
case 10 

worksheet = 'Miscellaneous'; 
outputfile = 'MiscellaneousOut.txt'; 
case 11 

worksheet = 'Oscillation'; 
outputfile = 'OscillationOut.txt'; 
end 

% Import the test filenames and test conditions from the selected 




% worksheet. 

[seriesnum,seriestext] = xlsread('MIT Test Plan',worksheet); 
numtests = size(seriesnum,l)-2; 

if series ==11 % For an unknown reason, the oscillation test worksheet 
% imports differently from the other worksheets, 
numtests =size(seriesnum,l); 
for testindex = 1 mumtests 

run(testindex) = (seriestext(testindex+2,3)); 

DandS(testindex,l) = seriesnum(testindex,l); 

DandS(testindex,2) = seriesnum(testindex,2); 
end 
else 

for testindex = 1 mumtests 
run(testindex) = (seriestext(testindex+2,3)); 

DandS(testindex,l) = seriesnum(testindex+2,l); 
DandS(testindex,2) = seriesnum(testindex+2,2); 
end 
end 

% Display the first and last files to be analyzed. 
firsttest=char(run(l)); 
lasttest=char(run(numtests)); 

fprintf('Autoanalyze will proceed from %s to %s \n', firsttest,lasttest) 

% Open the output file, 
warning off 
delete(outputfile) 
warning on 

manyrowsfid=fopen(outputfile,'a'); 

% Display which file is being processed and process the test file. Send 
% the test condition data that can not be extracted for the data file, 
for fileindex = 1 mumtests 
fhame = char(run(fileindex)); 
fprintf('Processing %s\n',fname) 

Depth = DandS(fileindex,l); 

Submergence = DandS(fileindex,2); 
if series == 8 % Steady force tests. 

Yaw Angle = seriesnum(fileindex+2,3); 

PitchAngle = seriesnum(fileindex+2,4); 

Rudder = 0; 

StemPlanes = 0; 
analyzemodxlssf; 
else % All unsteady motion tests. 




Yaw Angle = 0; 
PitchAngle = 0; 
Rudder = 0; 
StemPlanes = 0; 
analyzemodxls 
end 


end 

% Display that processing is complete and close the output file. 
fprintf('Processing Complete'); 
status = fclose(manyrowsfid); 




Appendix D. Sample Data File 

This is the first portion of the data from experiment OSC07. OSC07 involved 
oscillations in sway with forward velocity. The two rows of text at the top were added by 
the author for clarification and are not actually part of the data file. Not shown are the 
additional rows at the end of the data file which contain the ordered parameters for the 


experiment. 


1 Forces and Moments i 

Wave Height 
Channels (Not Used) 

[ Position Data in Controller Counts | 

1 Time Data I 

Alternate Time Data 

SSSfiSlCS 


IX |Y 


iRoli ! Pitch 1 

iTime Data! Dal 

te OatajOeltaT | 

Time DatelOettaT 

0.817 

- 0.031 

0.086 

0.086 

0.102 

- 0.115 

0.002 

0.001 

165 

316744 

798231 

503551 

13916 

37142.76 

37725 

0 

0.027 

0 

0.322 

- 0.099 

0.005 

0.021 

- 0.456 

- 0.026 

0 

0.001 

166 

319080 

798231 

503551 

13916 

37142.8 

37725 

0.039063 

0,073 

0.046 

0.667 

- 0.057 

0.005 

0.072 

- 0.796 

0.002 

- 0.001 

0.001 

227 

321424 

798231 

503551 

13916 

37142.84 

37725 

0.039063 

0.116 

0.043 

0.361 

- 0.006 

0.014 

- 0.091 

0.784 

- 0.079 

0.002 

0.001 

416 

323714 

798231 

503551 

13916 

37142.88 

37725 

0.042969 

0.156 

0.04 

- 0.052 

- 0.016 

0.005 

0.149 

- 0.492 

- 0.045 

0.001 

0.002 

636 

326137 

798230 

503551 

13916 

37142.92 

37725 

0.039063 

0.198 

0.042 

0.413 

- 0.036 

0 

0.584 

- 0,238 

0.012 

0.001 

0.003 

1350 

327779 

798231 

503551 

13916 

37142,96 

37725 

0.039063 

0,241 

0.043 

0.638 

- 0.074 

0.003 

0.347 

- 0.371 

- 0.079 

0.001 

0.001 

2548 

330343 

798230 

503551 

13916 

37143 

37725 

0.039063 

0.269 

0.028 

0.294 

- 0.086 

- 0.005 

0.489 

0.129 

- 0.111 

0 

0.001 

4159 

332879 

798232 

503551 

13916 

37143.04 

37725 

0.039063 

0.31 

0.041 

- 0.027 

- 0.046 

0.005 

0.295 

- 0.067 

- 0.049 

0.001 

0 

6398 

335523 

798230 

503551 

13916 

37143.08 

37725 

0.042969 

0.352 

0.042 

0.229 

- 0.101 

0.008 

- 1.175 

0.169 

- 0,065 

0.001 

0.001 

9262 

338054 

798231 

503551 

13916 

37143.12 

37725 

0.039063 

0.395 

0.043 

- 0.008 

0.011 

0.013 

- 1.092 

0.119 

- 0.004 

0 

0 

12747 

340486 

798230 

503551 

13916 

37143.16 

37725 

0.039063 

0.437 

0.042 

0.286 

0.024 

0.004 

0.519 

- 0.023 

0.033 

0 

0.001 

16652 

342728 

798231 

503551 

13916 

37143.2 

37725 

0.039063 

0.477 

0.04 

0.345 

- 0.059 

- 0.005 

0.794 

- 0.127 

- 0.118 

- 0.006 

0 

21057 

344921 

798230 

503551 

13916 

37143.24 

37725 

0.042969 

0.518 

0.041 

0.043 

- 0.167 

0.004 

0.646 

0.224 

- 0.166 

0 

0.001 

25744 

347005 

798231 

503551 

13916 

37143.28 

37725 

0.039063 

0.558 

0.04 

- 0.057 

0.042 

0.011 

0.361 

- 0.055 

0.07 

0 

0.002 

29136 

348342 

798231 

503551 

13916 

37143.32 

37725 

0.039063 

0.587 

0.029 

- 0.341 

- 0.067 

0.015 

- 0.454 

0.676 

0085 

0 

0.002 

34368 

350165 

798230 

503552 

13916 

37143.36 

37725 

0.042969 

0.629 

0.042 

- 0.234 

0.013 

0.011 

- 0.067 

0.159 

0026 

0.001 

0,007 

39725 

351745 

798232 

503551 

13916 

37143.4 

37725 

0.039063 

0.672 

0.043 

- 0.069 

0.184 

0.015 

0.009 

- 0.448 

- 0.033 

0 

0.002 

45221 

353150 

798231 

503551 

13916 

37143.44 

37725 

0.039063 

0.714 

0.042 

0.208 

- 0.013 

0.031 

- 0.441 

0.257 

- 0,029 

0.001 

0.001 

50493 

354288 

798231 

503551 

13916 

37143.48 

37725 

0.039063 

0.754 

0.04 

- 0.189 

- 0.152 

0.012 

0.088 

- 0.165 

0.004 

0.001 

0.001 

55941 

355252 

798232 

503551 

13916 

37143.52 

37725 

0.042969 

0.797 

0.043 

- 0.843 

- 0.053 

0.027 

0.165 

0.834 

0.084 

- 0.001 

0 

61351 

355897 

798232 

503551 

13916 

37143.56 

37725 

0.039063 

0.839 

0042 

- 0.409 

- 0.003 

0.017 

0.167 

- 0.171 

0-021 

0.001 

0 

64685 

356137 

798231 

503551 

13916 

37143.6 

37725 

0.039063 

0.88 

0.041 

- 0.109 

0.034 

0.014 

0.562 

0.097 

0.039 

0.001 

0.001 

70086 

3 K 319 

798231 

503551 

13916 

37143.64 

37725 

0.039063 

0.908 

0.020 

- 0.173 

- 0.03 

0.01 

- 0.12 

0.203 

- 0.003 

0.001 

0 

75293 

356303 

798231 

503551 

13916 

37143.68 

37725 

0.042969 

0.949 

0,041 

- 0.205 

- 0.018 

0.007 

- 0.403 

- 0.553 

0.049 

0 

0 

80797 

356089 

798231 

503551 

13916 

37143.72 

37725 

0.039063 

0.991 

0,042 

- 0.698 

- 0.089 

0.013 

0.284 

0.863 

0.061 

0.001 

0.001 

86001 

355606 

798232 

503551 

13916 

37143.76 

37725 

0.039063 

1.031 

0.04 

- 0.97 

0.019 

0.015 

0.106 

0.113 

- 0.007 

0.001 

0.001 

91428 

354722 

798232 

503551 

13916 

37143.8 

37725 

0.039063 

1,074 

0.043 

- 0.219 

- 0.052 

0.012 

0.192 

- 0.204 

0.035 

0 

0 

96569 

353570 

798231 

503551 

13916 

37143.84 

37725 

0.042969 

1.114 

0.04 

- 0.267 

0.033 

0.014 

0.26 

0.544 

0.066 

0 

0.002 

101985 

352184 

798231 

503551 

13916 

37143.88 

37725 

0.039063 

1.157 

0.043 

- 0.382 

0.05 

0.065 

- 0.104 

- 0.105 

- 0,024 

0 

0 

107466 

350626 

798231 

503551 

13916 

37143.92 

37725 

0,039063 

1.199 

0,042 

- 0.141 

- 0.035 

0.012 

- 0.241 

- 0.123 

0.036 

0.001 

0 

111116 

349500 

798231 

503551 

13916 

37143.96 

37725 

0.039063 

1.228 

0.029 

- 0.668 

0.002 

0.019 

0.342 

0.669 

0.02 

0001 

0 

116324 

347693 

798231 

503551 

13916 

37144 

37725 

0.042969 

1.268 

0.04 

- 0.688 

0.077 

0.014 

0.321 

- 0.271 

0.059 

0 

0.001 

121748 

345560 

798231 

503551 

13916 

37144.04 

37725 

0.039063 

1.311 

0.043 

- 0.136 

- 0.047 

- 0.004 

0.218 

0.261 

- 0.003 

0.001 

0.001 

127162 

343249 

798231 

503551 

13916 

37144.08 

37725 

0.039063 

1.353 

0.042 

- 0.19 

0.063 

0.003 

0.39 

- 0.264 

0.069 

0.001 

0.001 

132580 

340876 

798231 

503551 

13916 

37144.12 

37725 

0.039063 

1.395 

0.042 

- 0.062 

- 0.066 

0.008 

- 0.456 

- 0.179 

0.039 

0,001 

0.001 

137793 

338568 

798231 

503551 

13916 

37144.16 

37725 

0.042969 

1.436 

0,041 

- 0.208 

- 0.041 

0.011 

- 0.45 

- 0.086 

-0 018 

0 

0.001 

143271 

336076 

798231 

503551 

13916 

37144.2 

37725 

0.039063 

1.478 

0,042 

- 0.6 

- 0.094 

0.013 

0.581 

0.243 

0.072 

0,001 

0.001 

147161 

334225 

798232 

503551 

13916 

37144.24 

37725 

0.039063 

1.509 

0.031 

- 0.333 

0.107 

0.009 

- 0.048 

0.169 

-0 048 

■ 0 

0 

152338 

331709 

798232 

503551 

13916 

37144.28 

37725 

0.039063 

1.563 

0.054 

- 0.061 

- 0.154 

0.01 

- 0.065 

- 0.001 

- 0.025 

- 0.003 

0.001 

157760 

329060 

798231 

503551 

13916 

37144.32 

37725 

0.042969 

1.592 

0.029 

0.509 

0.079 

0.001 

0.307 

-0 365 

- 0.025 

0.002 

0 

163177 

326519 

798231 

503551 

13916 

37144.36 

37725 

0.039063 

1.634 

0,042 

0.11 

- 0.038 

0.001 

- 0.117 

-032 

0 117 

0 

0.002 

168907 

323971 

798231 

503551 

13916 

37144.4 

37725 

0.039063 

1.679 

0.045 

- 0.183 

- 0.019 

0.015 

- 0.508 

-0 109 

0.009 

0,001 

- 0.001 

174105 

321717 

798232 

503551 

13916 

37144.44 

37725 

0.039063 

1.719 

0.04 

- 0.264 

- 0.197 

0.014 

0.302 

0.189 

- 0.049 

0.001 

0.001 

179560 

319427 

798232 

503551 

13916 

37144.48 

37725 

0.042969 

1.761 

0042 

- 0.124 

- 0.011 

0.008 

0.177 

0.261 

- 0.152 

0.001 

0 

183179 

317973 

798232 

503551 

13916 

37144.52 

37725 

0.039063 

1.79 

0.029 

0.584 

- 0.259 

0.022 

0.136 

- 0.581 

0.012 

0 001 

0 

188340 

316043 

798231 

503551 

13916 

37144.56 

37725 

0.039063 

1.83 

0.04 

0.397 

0.078 

0.005 

0.792 

0.021 

0 092 

0 

0.001 

193773 

314203 

798231 

503551 

13916 

37144.6 

3772 S 

0.039063 

1.873 

0.043 

0.073 

- 0.127 

0.005 

- 0.214 

- 0.107 

-01 

0 

0.001 

198972 

312694 

798231 

503551 

13916 

37144.64 

37725 

0.042969 

1.913 

0.04 

0.153 

- 0.039 

0.013 

- 0.529 

0.065 

-012 

-0 001 

0.002 

204427 

311307 

798232 

503551 

13916 

37144.68 

37725 

0.039063 

1.956 

0.043 

0.091 

- 0.033 

0.012 

- 0.027 

- 0.299 

0073 

0.001 

0.001 

209887 

310112 

798232 

503551 

13916 

37144.72 

37725 

0.039063 

1.998 

0.042 

0.572 

- 0.008 

0.023 

- 1.172 

- 0.164 

- 0,008 

0,002 

0.001 

215065 

309219 

798232 

503551 

13916 

37144.76 

37725 

0.042969 

2.038 

0,04 

0.474 

- 0.261 

0.019 

- 0.152 

0.118 

- 0.047 

0.001 

0 

220485 

308554 

798232 

503551 

13916 

37144.8 

37725 

0.039063 

2.081 

0.043 

0.308 

0.067 

0.012 

0.431 

- 0.427 

0.049 

0.001 

0 

223845 

308301 

798232 

503551 

13916 

37144.84 

37725 

0.039063 

2.107 

0.026 

0.676 

- 0.112 

0.01 

0.179 

- 0.08 

- 0.019 

0.001 

0 

229289 

308112 

798231 

503551 

13916 

37144.88 

37725 

0.039063 

2.15 

0043 

0.12 

0.135 

0.017 

0.102 

0.201 

-0069 

0.001 

0 

234754 

308141 

798231 

503551 

13916 

37144.92 

37725 

0.042969 

2.192 

0.042 

0.345 

- 0.029 

0.014 

0.504 

- 0.113 

- 0.03 

0.002 

0.001 

239953 

308375 

798231 

503551 

13916 

37144.96 

37725 

0.039063 

2.233 

0.041 

0.592 

- 0.043 

0.011 

- 0.115 

- 0.098 

- 0.062 

0.001 

0.001 

245392 

308866 

798231 

503551 

13916 

37145 

37725 

0.039063 

2,275 

0.042 

0.997 

- 0.122 

0.026 

0.062 

- 0.601 

- 0.046 

- 0.001 

0 

250823 

309761 

798231 

503551 

13916 

37145.04 

37725 

0.039063 

2.317 

0.042 

0.852 

0.007 

0.019 

0.266 

0.224 

-0 083 

- 0.001 

0.001 

256250 

310980 

798231 

503551 

13916 

37145.08 

37725 

0.039063 

2.36 

0.043 

0.216 

- 0.065 

0.013 

- 0 . 16 B 

- 0.435 

0.005 

0.001 

0.002 

261434 

312375 

798231 

503551 

13916 

37145.12 

37725 

0.042969 

2.4 

0.04 

0.158 

- 0.005 

0.02 

- 0.38 

0.164 

- 0.049 

- 0.021 

0.003 

265073 

313395 

798231 

503551 

13916 

37145.16 

37725 

0.039063 

2.429 

0029 

0.509 

- 0.084 

0.013 

0.315 

• 0.457 

- 0.018 

- 0.001 

0.001 

270532 

315070 

798231 

503551 

13916 

37145.2 

37725 

0.039063 

2,483 

0.054 

0.697 

- 0.043 

0.016 

0.169 

0.255 

-0115 

0.001 

0.003 

275725 

316843 

798232 

503551 

13916 

37145.24 

37725 

0.039063 

2.512 

0.029 

0.476 

- 0.116 

0.016 

- 0.246 

- 0.355 

- 0.014 

0 

V 0 

280904 

318874 

798232 

503551 

13916 

37145.28 

37725 

0.042969 

2.552 

0.04 

0.748 

- 0.047 

0.021 

- 0.064 

- 0.639 

0.021 

0.001 

' 0 

286336 

321203 

798232 

503551 

13916 

37145.32 

37725 

0.039063 

2.594 

0.042 

0.305 

- 0.001 

0.051 

- 0.252 

0,825 

-0 066 

0 

- 0.001 

291770 

323613 

798232 

503551 

13916 

37145.36 

37725 

0.039063 

2.637 

0.043 

- 0.091 

0.014 

0.012 

- 0.507 

-0 449 

• 0.05 

0.001 

0.001 

297219 

326056 

798232 

503551 

13916 

37145.4 

37725 

0.042969 

2.679 

0.042 

0.382 

- 0.016 

0.021 

0.265 

- 0.11 

- 0.031 

0.001 

0.001 

300600 

327567 

798232 

503551 

13916 

37145.44 

37725 

0.039063 

2.72 

0.041 

0.648 

- 0.076 

0.012 

0.397 

- 0.193 

- 0.112 

0.001 

0 

306050 

330104 

798232 

603551 

13916 

37145.48 

37725 

0.039063 

2.762 

0.042 

0.426 

- 0.126 

0.025 

0.019 

- 0,047 

- 0.095 

0 

0.001 

311488 

332753 

798232 

503551 

13916 

37145.52 

37725 

0.039063 

2.791 

0.029 

0.012 

- 0.131 

0.016 

0.711 

-0 134 

- 0.025 

0 

0.001 

316927 

335402 

798232 

503550 

13916 

37145.56 

37725 

0.042969 

2.833 

0.042 

0.175 

0.014 

0.007 

0.311 

0.11 

- 0.062 

0.001 

0.001 

322098 

337845 

798232 

503551 

13916 

37145.6 

37725 

0.039063 

2.874 

0,041 

- 0.033 

- 0.119 

- 0.018 

- 0.181 

0.614 

• 0.076 

0.001 

0.001 

327555 

340282 

798232 

503551 

13916 

37145.64 

37725 

0.039063 

2916 

0.042 







Appendix E. AnalyzemodXlsSF.m 

% AnalyzemodXlsSF'.m 
% Erik Oiler, 2003 

% Performs analyses for Steady Force Tests only 

% This program is designed to be called by AutoAnalyzeXls.m and requires 
%no manual intervention. 

% Initialize matrices. 

[Data] = 0; 

[Datain] =0; 

% Get the input data and find the length of the file, 
infname = strcat(fname,'.xls'); 

existencecheck = exist(infname); % Determines if the input file exists, 
if existencecheck ==0 

fprintf('%s does not exist.\n',infname) 
return 
end 

Datain = xlsread(infiiame); % Reads the input file. 

NumLines = size(Datain,l)-4; % Gets the number of data samples in 
% the input file. 

% Gets the date and converts from EXCEL to MATLAB format. 
Date=datestr(Datain(1,15)-36525,2); 

% Build a date-time string to ensure the correct gains are applied. 

% Add 693960 to the date to get number of days from 0000. 

% Divide the number of minutes by 86400 to get fractional days. 

% MATLAB date serial numbers are in the form; 

% days since OOOO.fraction of a day 
DateNum=Datain(l, 15)+693960+Datain(l, 14)/86400; 

% Deterimine the ordered velocity. 

XSpd = Datain(NumLines+3,10); 
if isnan(XSpd)=l 
XSpd = 0; 
end 

% Determine the ordered distance of travel. 

XDist = Datain(NumLines+3,l 1); 
if isnan(XDist)==l 



XDist = 0; 
end 

% Determine the ordered length of time of the experiment. 

Time = Datain(NumLines+4,5); 

% Determine actual sample frequency based upon the measured interval 
% between data points. Necessary because the control software did not 
% always sample at the ordered sample rate. 

SampleFreq = l/mean(Datain(l:NumLines,16)); 


% Drop the first 1.2 seconds of data to account for acceleration. At 

# % the sample frequency of 25 Hz, drop the first 30 points. Shift the 
% other points up in the array. 

dropgap = round(l .2*SampleFreq); 
for I = 1 :NumLines-dropgap; 

Datain(I,l: 16)=Datain(I+dropgap, 1:16); 

# end 

NumLines = NumLines - dropgap; 

% Drop data recorded after the model stopped moving along the X-axis 
^ % for non-Inertial tests. 

NumLines2=NumLines; 
for 1 = NumLines-1 :-l :75 

if (Datain(I,9)==Datain(I+l,9)) %Looks for constant X position. 

NumLines2 = I; 
end 

® end 

NumLines = NumLines2; 

% Determine the time of the first data point. Used for converting 
% from time past midnight to time of run. 

# TimeStart = Datain( 1,14); 


% The inverse sensitivity matrix (B). The sensitivity matrix used 
% depends on the when the experiment was performed. Ealry experiments 
% used a different load cell than later ones. 

if((datenum(Date)>=datenum(2002,07,17))&... 
(datenum(Date)<datenum(2002,10,31))) 

[B] = [ 0.3901 0.0029 -0.0071 -0.0010 -0.0024 -0.0016; 

0.0023 0.3887 0.0016 0.0025-0.0039 0.0019; 

0.0139 0.0129 1.5006-0.0002-0.0215 -0.0018; 

0.0000-0.0001 0.0016 0.0069 0.0000-0.0001; 
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-0.0001 0.0000 0.0018 0.0000 0.0070 0.0000; 

-0.0003 -0.0002 0.0000 0.0000 0.0000 0.0108]; 
else %valid on and after Oct 31, 2002 

[B] = [ 0.3856 0.0020 -0.0016 -0.0008 0.0001 -0.0030; 

0.0023 0.3811 -0.0040 0.0012-0.0034 0.0031; 

0.0093 0.0024 1.5109 0.0068-0.0231 -0.0014; 

0.0001 0.0000 0.0013 0.0069 0.0000 0.0000; 

0.0000 0.0000 0.0007 0.0000 0.0069-0.0001; 

-0.0003 0.0001 0.0000 0.0000 0.0000 0.0106]; 

end 

% Establish Gains and Excitation Voltage 
if (DateNum>731434.61458) % August 6, 2002 1445. 

[Gain] = [4000;4000; 1000;! 000; 1000;4000]; % y z x pitch yaw roll 
[Vexc] = [10; 10; 2.5; 10; 10; 10]; 
end 

% Calculate the Conversion Factors (CF) between voltage and force 
% or moment 

CFtemp = Gain.*Vexc*10'^-6; 

CF = zeros(6,6); 

CF(1,1)= CFtemp(l); 

CF(2,2)- CFtemp(2); 

CF(3,3)= CFtemp(3); 

CF(4,4)= CFtemp(4); 

CF(5,5)= CFtemp(5); 

CF(6,6)= CFtemp(6); 

% Establish multipliers to convert from controller counts to MKS units. 
% Multipliers are different for different time intervals due to system 
% upgrades. 

% Between 0800 Oct 24,2002 and 0800 Dec 1, 2002. 
if ((DateNum>731513.33333)&(DateNum<731551.33333)) 
xfactor = 3850; 
yfactor = 2410; 
zfactor = 2114; 
pitchfactor = 972; 
yawfactor= 1818; 

elseif DateNum>731551.33333 % After 0800 Dec 1,2002 
xfactor = 3850; 
yfactor = 2410; 
zfactor = 2114; 
pitchfactor =155; 
yawfactor= 1818; 
else % Before Oct 24, 2003 





xfactor = 3850; 
yfactor = 2410; 
zfactor = 4921; 
pitchfactor =1111; 
y awfactor = 1025; 
end 

% Convert from voltages to forces 
for I = 1 : NumLines 

Datain(I,l:6) = (CF^-l*B*Datain(I,l;6)')'; 

Datain(I,9) =Datain(I,9) /xfactor; 

Datain(I,10) = Datain(I,10) / yfactor; 

Datain(I,l 1) = Datain(I,l 1) / zfactor; 

Datain(I,12) = Datain(1,12) /y awfactor; 

Datain(I,13) = Datain(I,13) / pitchfactor; 

DatainP,14) = Datain0,14) - TimeStart; 
end 

% Find the mean force 
forJ=l:6 

MeanForce(J) = real(mean(Datain(l :NumLines,J))); 
end 

% PROCESSING STAGE 

% Shift the origin of the coordinate system from the origin of the load 
% cell to vessel amidships, 
shiftlength = 0.0522; % meters 

MeanForce(4) = MeanForce(4) - MeanForce(2) * shiftlength; % Pitch 
MeanForce(5) = MeanForce(5) + MeanForce(l) * shiftlength; % Yaw 

% Print the results to a common row output file. 
fprintf(manyrowsfid,'%s \t %s\t %6.2f,ftiame,Date,Depth); 
^rintf(manyrowsfid,'\t%6.2f \t%5.1 f \t %5.2f\t%5 .If,... 

Submergence,XSpd,XDist,Time); 
fyrintf(manyrowsfid,' \t %2.0f \t%2.0f\t %3.1f\t %3.1f,... 

PitchAngle,Y awAngle,StemPlanes,Rudder); 
fprintf(manyrowsfid,' \t %7.4f,... 

MeanForce(3),MeanForce(l),MeanForce(2),MeanForce(6),MeanForce(4),. 

MeanForce(5)); 

fprintf(manyrowsfid,'\n'); 



Appendix F. Output File from AnalyzemodXlsSF.m called 

• by AutoanalyzeXls.m for Steady Force Tests 

This is part of the output file from AnalyzemodXlsSF.m when called by 
AutoanalyzeXls.m to analyze steady force tests. The two rows of text at the top were 
added by the author for clarification and are not actually part of the data file. When using 
this table please note the following: 

# SF67-72, 80-90, 144-185 were performed with control surfaces removed. 

SF73 was a test to determine the effect of system noise with no motion. 

SF74 was performed with no motion but with a 4 oz weight attached near midships. 

SF75 was performed with no motion but with a 2 oz weight attached near midships. 



Date 


lESSSHi 

Spd 

Xdist 

Time 

Pitch 

Yaw 

SP 


X 

Y 

Z 

Ron 

Pitch 

Yaw 

1 


cm 

mm 

cm/s 

cm 

s 



Degrees 

Degrees 

N 

N 

N 

N-m 

N-m 

N-m 



SF6t 

03/05/03 

79.5 

356.2 

34 

700 

22 

■El 


0 

0 

6.7506 

-0.5651 

0.736 

-0.001 

-0.0746 

-0.0205 

SF62 

omsioz 

79.5 

356.2 

34 

700 

22 

0 

4 

0 

0 

10.7151 

-0.1987 

0.3827 

0.0051 

-0.032 

0.0089 

SF63 

03/05/03 

79.5 

3562 

34 

700 

22 

0 

8 

0 

0 

10.9515 

-0.1364 

0.4721 

0.0067 

-0.0244 

0.0176 

SF64 

03/05/03 

79.5 

356,2 

68 

700 

11 

0 

0 

0 

0 

5.3257 

-0.5732 

0.7378 

-0.0013 

-0.0782 

-0.0175 

SF65 

03/05/03 

79.5 

356.2 

68 

700 

11 

0 

4 

0 

0 

7.3517 

0.0361 

0.7981 

-0.0072 

-0.0851 

0.0612 

SF66 

03/05/03 

79.5 

356.2 

68 

700 

11 

0 

8 

0 

0 

6.7879 

0.3516 

0.6998 

-0.0046 

-0.0372 

0,1184 

SF67 

03/05/03 

79.5 

356.2 

34 

700 

22 

0 

0 

0 

0 

19.8251 

0.6296 

-0.2936 

-0.0133 

-0.0042 

0.0361 

SF66 

03/05/03 

79.5 

356.2 

34 

700 

22 

0 

4 

0 

0 

20.801 

0.3637 

-0.0698 

-0,0069 

-0.0256 

0.0391 

SF69 

03/05/03 

79.5 

356.2 

34 

700 

22 

0 

8 

0 

0 

17.9008 

0.6628 

-0.175 

-0.0165 

-0.0003 

0,0694 

SF70 

03/05/03 

79.5 

356.2 

68 

700 

11 

0 

0 

0 

0 

18.9798 

0.5325 

0.1027 

-0.0176 

-0.0432 

0.0359 

SF71 

03^)5/03 

79.5 

356.2 

68 

700 

11 

0 

4 

0 

0 

17.655 

0.5197 

-0.1675 

-0.0142 

-0.0002 

0.1008 

SF72 

03^5/03 

79.5 

356.2 

68 

700 

11 

0 

8 

0 

0 

17.5917 

0.9081 

-0.0402 

-0.017 

0.0217 

0.1842 

SF73 

03/05/03 

79.5 

356.2 

0 

0 

20 

0 

0 

0 

0 

-3.8807 

0.1229 

-0.1337 

-0.0056 

0.0123 

0.0093 

SF74 

03/07/03 

79.5 

356.2 

0 

0 

20 

0 

0 

0 

0 

17.5033 

-0.1331 

1.0175 

0.0292 

-0.0721 

-0.014 

SF75 

03/07/03 

79.5 

356,2 

0 

0 

20 

0 

0 

0 

0 

18.2672 

-0.3689 

0.961 

0,0215 

-0.054 

-0.0237 

SF76 

04/04/03 

79.5 

356,2 

100 

700 

10 

0 

0 

0 

0 

14.4912 

-0.0623 

0.3053 

0.0178 

-0.098 

-0.0183 

SF77 

04/04/03 

79.5 

3562 

100 

700 

10 

0 

4 

0 

0 

7.3252 

0.2929 

0.257 

0.0059 

-0.0875 

00822 

SF7e 

04/04/03 

79.5 

356.2 

100 

700 

10 

0 

8 

0 

0 

7.8204 

0.5732 

0,469 

0.0035 

-00797 

01728 

SF78B 

04/04/03 

79.5 

356.2 

100 

700 

10 

0 

12 

0 

0 

6.3188 

1.2356 

0.8194 

0.0007 

-0.0858 

0.2656 

SF79 

04/04/03 

79.5 

356.2 

100 

700 

10 

0 

16 

0 

0 

5.8774 

2.1367 

0.0009 

-0.0019 

-0.1104 

0.3389 

SF80 

04/04/03 

79.5 

356.2 

100 

700 

10 

0 

0 

0 

0 

8.6834 

0.0818 

0.5544 

-0.0087 

-0.0791 

-0.0039 

SF81 

04/04/03 

79.5 

356.2 

100 

700 

10 

0 

4 

0 

0 

12.1557 

0.4921 

0.4522 

-0.0008 

-0.0802 

0.1348 

SF82 

04/04/03 

79 5 

356.2 

100 

700 

10 

0 

8 

0 

0 

7.9449 

1.0495 

0.0418 

-0.0179 

-0 0306 

0.2852 

SF82B 

04/04/03 

79.5 

356.2 

100 

700 

10 

0 

12 

0 

0 

10.9136 

1.0141 

0.4823 

0 

-0,1804 

03774 

SF83 

04/04/03 

79.5 

356.2 

100 

700 

10 

0 

16 

0 

0 

8.0716 

2.2319 

-0.2717 

-0.0166 

-0.0772 

0.5317 

SF84 

04/04/03 

79.5 

356.2 

100 

700 

10 

0 

0 

0 

0 

23.3779 

-0.1231 

0.3215 

0.0153 

-0.0857 

-0 0266 

SF85 

04/04/03 

79.5 

356.2 

100 

700 

10 

0 

4 

0 

0 

26.1843 

0.1264 

0.1717 

0.016 

-0.0563 

0.1104 

SF86 

04/04/03 

79.5 

356 2 

100 

700 

10 

0 

8 

0 

0 

22.1398 

0.4683 

0.5599 

0.0112 

-0.0638 

0.2422 

SF87 

04/04/03 

79.5 

356 2 

100 

700 

10 

0 

12 

0 

0 

20.9844 

0.9228 

0.5964 

0.0115 

-0 0872 

0,3704 

SF88 

04/04/03 

79.5 

356,2 

100 

700 

10 

4 

0 

0 

0 

21.7807 

-0.3638 

0.1393 

0.0182 

-0 0021 

-0.0355 

SF89 

04/04/03 

795 

^ 356 2 

100 

700 

10 

8 

0 

0 

0 

23.6841 

-0.569 

-0.1426 

0,0206 

0 1086 

-0.0427 

SF90 

04/04/03 

79.5 

356 2 

100 

700 

10 

12 

0 

0 

0 

21.4243 

-0.475 

-0.5479 

0.0202 

0.2067 

-0 036 

SF91 

04/04/03 

79.5 

356 2 

100 

700 

10 

0 

0 

0 

0 

7.5073 

0.0728 

-0.015 

-0,0008 

-0,1033 

-0.0076 

SF92 

04/04/03 

795 

356,2 

100 

700 

10 

0 

4 

0 

0 

10.8973 

0.7809 

-0.1784 

-0.0129 

-0.0639 

0,1026 

SF93 

04/04/03 

79.5 

356.2 

10O 

700 

10 

0 

8 

0 

0 

6.0916 

1.2924 

0.3105 

-0.0203 

-0.0868 

0.2104 

SF94 

04/04/03 

79.5 

356.2 

100 

700 

10 

0 

12 

0 

0 

8.3035 

1.6519 

0.2402 

-0.0108 

-0.078 

02797 

SF95 

04/04/03 

79.5 

356.2 

100 

700 

10 

4 

0 

0 

0 

6.013 

-0.0121 

0.0331 

-0,0065 

-0.0191 

-0.0095 

SF96 

04/04/03 

79.5 

356.2 

100 

700 

10 

8 

0 

0 

0 

8.3074 

0.0053 

-0.3269 

-0,0063 

0.0113 

-0.0062 

SF97 

04/04/03 

79.5 

356.2 

100 

700 

10 

12 

0 

0 

0 

7.1199 

0.0485 

-0.993 

-0.0044 

0.0543 

-0.0039 

SF98 

04/07/03 

79.5 

356.2 

100 

700 

10 

0 

0 

0 

0 

14.3126 

-0.6259 

0.5736 

0,0084 

-0.1263 

-0.0457 

SF99 

04/07/03 

79.5 

356,2 

75 

700 

13 

p 

0 

0 

0 

17.1046 

-0.6789 

0.2246 

0.0136 

-0,0526 

-0.0474 

SF100 

04/07/03 

79.5 

210.9 

100 

700 

10 

b 

0 

0 

0 

28.6384 

-0.3125 

-0.3894 

0.0441 

-0 0554 

-0 0408 

SF101 

04/07/03 

79.5 

210.9 

75 

700 

13 

0 

0 

0 

0 

21.3003 

-0.3711 

-0.4691 

0.0403 

-0,0416 

-0.0368 

SF102 

04/07/03 

79.5 

501.5 

100 

700 

10 

0 

0 

0 

0 

8.4768 

0.2046 

0.2218 

-0.0085 

-0.1139 

0 

SF103 

04/07/03 

79.5 

501 5 

75 

700 

13 

0 

0 

0 

0 

8.104 

0.2389 

-0.1237 

-0.0059 

-0.0566 

0.0065 

SF104 

04/07/03 

79.5 

356.2 

100 

700 

10 

0 

4 

0 

0 

0.0083 

-0.1356 

0.7407 

-0.0239 

-0 1217 

0.0659 

SF105 

04/07/03 

79.5 

356.2 

75 

700 

13 

0 

4 

0 

0 

-5.3657 

0.1665 

0.4406 

-0.0349 

-0.0367 

0.0572 

SF106 

04/07/03 

79.5 

210.9 

100 

700 

10 

0 

4 

0 

0 

7.9829 

0.5966 

-0.2151 

-0.0043 

-0.0736 

0.0959 

SF107 

04/07/03 

79.5 

210.9 

75 

700 

13 

0 

4 

0 

0 

11.5173 

0.4629 

-0.9813 

0.0007 

-0.0326 

0.0647 

SF108 

04/07/03 

79.5 

501.5 

100 

700 

10 

0 

4 

0 

0 

13.1946 

0.6244 

-0.1566 

0.0198 

-0.0778 

0.0915 

SF109 

04/07/03 

79.5 

501.5 

75 

700 

13 

0 

4 

0 

0 

17.8258 

0.6527 

-0.3395 

0.0142 

-0.0294 

0.0722 

SF110 

04/07/03 

79.5 

3562 

100 

700 

10 

0 

8 ' 

0 

0 

-7.1228 

0.7253 

0.6748 

-0.0415 

-0.1104 

0.1874 

SF111 

04/07/03 

79.5 

356 2 

75 

700 

13 

0 

8 

0 

0 

-4.6664 

0.2552 

0.8404 

-0.0352 

-0.044 

0.1106 

SF112 

04/07/03 

79.5 

210.9 

100 

700 

10 

0 

8 

0 

0 

5.9871 

1.1071 

0.4481 

-0.0174 

-0,0857 

0.199 

SF113 

04/07/03 

79.5 

210.9 

75 

700 

13 

0 

8 

0 

0 

9.4164 

0.8015 

0.934 

-0.0123 

-0.033 

0.1305 

SF114 

04/07/03 

79.5 

501.5 

100 

700 

10 

0 

8 

0 

0 

-5.3702 

0.4664 

0.8559 

-0.0415 

-0.1143 

0.1744 

SF115 

04/07/03 

79.5 

501.5 

75 

700 

13 

0 

8 

0 

0 

-5.1133 

0.0701 

0.808 

-0.0382 

-0.0898 

0.1038 

SF116 

04/07/03 

79.5 

356.2 

100 

700 

10 

0 

12 

0 

0 

-4,4508 

1.5926 

0.1168 

-0.0049 

-0.0714 

0.2824 

SF117 

04/07/03 

79.5 

356.2 

75 

700 

13 

0 

12 

0 

0 

-1.2336 

0.8619 

0.0513 

0.0016 

-0.0279 

0.1713 

SF118 

04/07/03 

79.5 

210.9 

100 

700 

10 

0 

12 

0 

0 

7.9814 

1.8835 

-0.1084 

0.0196 

-0.0386 

0.2878 

SF119 

04/07/03 

79.5 

210.9 

75 

700 

13 

0 

12 

0 

0 

9.5864 

1.4279 

0.1144 

0.0221 

0.0022 

0.1957 

SF120 

04/07/03 

79.5 

501.5 

100 

700 

10 

0 

12 

0 

0 

-7.2312 

1.2275 

0.894 

-0.0446 

-0.1012 

0.2688 

SF121 

04/07/03 

79.5 

501.5 

75 

700 

13 

0 

12 

0 

0 

1.4662 

0.8562 

0.0813 

0.0017 

-0.0339 

0.1749 

SF122 

04/07/03 

79.5 

356,2 

100 

700 

10 

4 

0 

0 

0 

1.5423 

-0.272 

-0.0287 

0.0086 

-0.0199 

-0.0161 

SF123 

04/07/03 

79.5 

356 2 

75 

700 

13 

4 

0 

0 

0 

-1.8708 

-0,1491 

0.0406 

0.0037 

-0.003 

-0.0077 

SF124 

04/07/03 

79.5 

210,9 

100 

700 

10 

4 

0 

0 

0 

11.5445 

0.3353 

-1.0374 

0.0309 

0.0335 

0.0038 

SF125 

04/07/03 

79.5 

210.9 

75 

700 

13 

4 

0 

0 

0 

14.0924 

0.2787 

-0.1434 

0.0275 

0.0337 

0,0086 








73 







Depth 

Subm. 

Spd 

Xdist 

Time 

1 Pitch 1 


HuaHi 


X 

Y 

Z 

Roll 





cm 

mm 

cm/s 

cm 

s 


1 Degrees 



N 

N 

N 

N-m 

IHZSjSHI 

HCSSHi 

SF126 

04/07/03 

79,5 

501.5 

100 

700 

10 

4 

0 

0 

0 

1.5582 

-0.3684 

-0.1921 

0.0073 

0.0015 

-0.0231 

SF127 

04/07/03 

79.5 

501.5 

75 

700 

13 

4 

0 

0 

0 

1.7731 

-0.0883 

-0.192 

0.0027 

0.0111 

-0.0039 

SF128 

04/07/03 

79.5 

356.2 

100 

700 

10 

8 

0 

0 

0 

-4.1236 

-0.1895 

^.6621 

0.0044 

0.0651 

-0.0037 

SF129 

04/07/03 

79 5 

356.2 

75 

700 

13 

8 

0 

0 

0 

-2.5511 

-0.2603 

-0.4944 

0.0051 

0.0613 

-0,0105 

SF130 

04/07/03 

79.5 

210.9 

100 

700 

10 

8 

0 

0 

0 

5.4903 

0.0547 

-0.8964 

0.0277 

0.0674 

0.0062 

SF131 

04/07/03 

79.5 

210.9 

75 

700 

13 

8 

0 

0 

0 

8.7811 

0.239 

-0.7814 

0,026 

0.0676 

0.0145 

SF132 

04/07/03 

79.5 

501.5 

100 

700 

10 

8 

0 

0 

0 

-1,9805 

-0.1559 

-0.6271 

-0.0014 

0.0606 

-0.0079 

SF133 

04/07/03 

79,5 

501.5 

75 

700 

13 

8 

0 

0 

0 

1.3509 

-0.0493 

-0.4612 

0.0006 

0.0524 

-0.0012 

SF134 

04/07/03 

79 5 

356.2 

100 

700 

10 

12 

0 

0 

0 

2.9365 

-0.4546 

-1.0079 

0.0183 

0.1172 

-0.0088 

SF135 

04/07/03 

79.5 

356.2 

75 

700 

13 

12 

0 

0 

0 

0.1206 

-0.2504 

-0.7061 

0.0112 

0.0867 

-0.0135 

SF136 

04/07/03 

79.5 

210.9 

100 

700 

10 

12 

0 

0 

0 

10.111 

0.2624 

-1.6039 

0.027 

0.1288 

0.0077 

SF137 

04/07/03 

79.5 

210.9 

75 

700 

13 

12 

0 

0 

0 

13.6603 

0.2164 

-1.1318 

0.0324 

0.1141 

0.0176 

SF138 

04/07/03 

79.5 

501.5 

100 

700 

10 

12 

0 

0 

0 

1.0817 

-0.526 

-0.8771 

0.0042 

0.1005 

-0.0276 

SF139 

04/07/03 

79.5 

501.5 

75 

700 

13 

12 

0 

0 

0 

0.7022 

-0.1462 

-0.7014 

0.0016 

0.0688 

-0.004 

SF140 

04/07/03 

79.5 

501.5 

100 

700 

10 

0 

8 

0 

0 

8.2723 

0.8434 

0.7545 

-0.0361 

-0.0988 

0.194 

SF141 

04/07/03 

79.5 

501.5 

75 

700 

13 

0 

8 

0 

0 

6.0705 

0,7525 

0.0117 

-0.006 

-0.0367 

0.134 

SF142 

04/07/03 

79,5 

210.9 

100 

700 

10 

8 

0 

0 

0 

14.5032 

0.172 

-0.7051 

0.0024 

0.0661 

0.0052 

SF143 

04/07/03 

79.5 

210.9 

75 

700 

13 

4 

0 

0 

0 

19.6126 

0.2523 

-0.4773 

0.005 

0.013 

0.0121 

SF144 

04/07/03 

79.5 

356.2 

100 

700 

10 

0 

0 

0 

0 

11.9977 

-0.3885 

0.1854 

0.0021 

-0.0888 

-0.0355 

SF145 

04/07/03 

79.5 

356 2 

75 

700 

13 

0 

0 

0 

0 

12.6364 

0.0928 

-0.0909 

-0.0048 

-0.0267 

-0.0034 

SF146 

04/07/03 

79.5 

210.9 

100 

700 

10 

0 

0 

0 

0 

11.2486 

0,0564 

-0.3442 

0.0016 

-0.0728 

-0.0127 

SF147 

04/07/03 

79.5 

210.9 

75 

700 

13 

0 

0 

0 

0 

12.6827 

0.0567 

-0.2317 

0.004 

-0.0304 

-0.008 

SF148 

04/07/03 

79.5 

501.5 

100 

700 

10 

0 

0 

0 

0 

7.8708 

-0.3903 

0.2886 

-0.0143 

-0.0798 

-0.0368 

SF149 

04/07/03 

79.5 

501.5 

75 

700 

13 

0 

0 

0 

0 

7.3999 

-0.2902 

0.1767 

-0.0129 

-0.0526 

-0.0216 

SF150 

04/07/03 

79.5 

356.2 

100 

700 

10 

0 

4 

0 

0 

-0.8426 

0.0721 

0.4963 

-0.0322 

-0.0894 

0.1153 

SF151 

04/07/03 

79.5 

356.2 

75 

700 

13 

0 

4 

0 

0 

2.0863 

-0.1088 

0.6661 

-0.0248 

-0.0585 

0.0663 

SF152 

04/07/03 

79.5 

210.9 

100 

700 

10 

0 

4 

0 

0 

8.6837 

0.2023 

0.3061 

-0.0066 

-0.082 

0.1125 

SF153 

04/07/03 

79.5 

210.9 

75 

700 

13 

0 

4 

0 

0 

11.8808 

0.2309 

-0.4769 

-0.0046 

-0.003 

0.0779 

SF154 

04/07/03 

79.5 

501.5 

100 

700 

10 

0 

4 

0 

0 

4.0575 

0.1175 

0.3141 

-0.0248 

-0.0839 

0.1181 

SF155 

04/07/03 

79.5 

501.5 

75 

700 

13 

0 

4 

0 

0 

11.2072 

0.0949 

0.188 

-0.0183 

-0.0416 

0.0777 

SF156 

04/07/03 

795 

356.2 

100 

700 

10 

0 

8 

0 

0 

1.8573 

0.3133 

0.6259 

-0.0283 

-0.0693 

0.2394 

SF157 

04/07/03 

79.5 

356.2 

75 

700 

13 

0 

8 

0 

0 

10.5516 

0.6441 

-0.2885 

0.0098 

0.0026 

0.1728 

SF158 

04/07/03 

79.5 

210.9 

100 

700 

10 

0 

8 

0 

0 

12.9229 

1.0197 

0.0742 

-0.0048 

-0.0546 

0.2681 

SF159 

04/07/03 

79.5 

210.9 

75 

700 

13 

0 

8 

0 

0 

11.4304 

0.4812 

0.8981 

-0.0097 

-0.0339 

0.1663 

SF160 

04/07/03 

79.5 

501.5 

100 

700 

10 

0 

8 

0 

0 

-2.0303 

0.4854 

0.2242 

-0.0035 

-0.0645 

0.253 

SF161 

04/07/03 

79.5 

501.5 

75 

700 

13 

0 

8 

0 

0 

4.6134 

0.7048 

0.1588 

-0.0186 

-0.0291 

0.1856 

SF162 

04/07/03 

79.5 

356,2 

100 

700 

10 

0 

12 

0 

0 

-5.1568 

1.2942 

0.4418 

-0.0165 

-0.0632 

0.4016 

SF163 

04/07/03 

79.5 

356.2 

75 

700 

13 

0 

12 

0 

0 

3.8291 

0.9787 

-0.1571 

0.0085 

0.0067 

0.2613 

SF164 

04/07/03 

79.5 

210.9 

100 

700 

10 

0 

12 

0 

0 

10.272 

1.5353 

0.1432 

-0.0012 

-0.038 

0.4029 

SF165 

04/07/03 

79.5 

210.9 

75 

700 

13 

0 

12 

0 

0 

13.0686 

1.065 

-0.8274 

0.0037 

-0.0076 

0.2577 

SF166 

04/07/03 

79.5 

501.5 

100 

700 

10 

0 

12 

0 

0 

-7.5024 

1.4947 

0.18 

-0.017 

-0.0564 

0.4133 

SF167 

04/07/03 

79.5 

501.5 

75 

700 

13 

0 

12 

0 

0 

6.2886 

0.9665 

-0.0814 

0.0012 

-0.0278 

0.2643 

SF168 

04/07/03 

79.5 

356.2 

100 

700 

10 

4 

0 

0 

0 

7.7978 

-0.6456 

0.5741 

-0.0051 

0.0123 

-0.0521 

SF169 

04/07/03 

79.5 

356.2 

75 

700 

13 

4 

0 

0 

0 

9.444 

-0.1815 

0.1395 

-0.0051 

0.0376 

-0.0216 

SF170 

04/07/03 

79.5 

210.9 

100 

700 

10 

4 

0 

0 

0 

3.7602 

0.8067 

-0.3934 

-0.025 

0.0252 

0.0285 

SF171 

04/07/03 

79.5 

210.9 

75 

700 

13 

4 

0 

0 

0 

3.9603 

0.7468 

0.0696 

-0.0282 

0.0341 

0.04 

SF172 

04/07/03 

79.5 

501.5 

100 

700 

10 

4 

0 

0 

0 

-3.5378 

-0.1568 

-0.0781 

0.0035 

0.0379 

-0.027 

SF173 

04/07/03 

79.5 

501.5 

75 

700 

13 

4 

0 

0 

0 

-3.392 

-0.1565 

0.047 

0.0039 

0.0278 

-0.0167 

SF174 

04/07/03 

79.5 

356.2 

100 

700 

10 

8 

0 

0 

0 

8.2774 

0.2489 

-0.6666 

0.009 

0.1331 

-0.0013 

SF175 

04/07/03 

79.5 

356.2 

75 

700 

13 

8 

0 

0 

0 

15.9178 

0.3335 

-0.6135 

0.0198 

0.1163 

0.0029 

SF176 

• 04/07/03 

79,5 

210.9 

100 

700 

10 

8 

0 

0 

0 

14.9037 

0.3619 

-0.6767 

0.0128 

0.1374 

0.0167 

SF177 

04/07/03 

79.5 

210.9 

75 

700 

13 

8 

0 

0 

0 

19.9008 

0.7155 

-0.5469 

0.0154 

0.1086 

0.0298 

SF178 

04/07/03 

79.5 

501.5 

100 

700 

10 

8 

0 

0 

0 

13.905 

-0.3111 

-0.2381 

0.0055 

0.1145 

-0.0191 

SF179 

04/07/03 

79.5 

501.5 

75 

700 

13 

8 

0 

0 

0 

12.3352 

-0.1277 

-0.1732 

0.002 

0.0754 

-0.0116 

SF180 

04/07/03 

79.5 

356 2 

100 

700 

10 

12 

0 

0 

0 

18.163 

0.2045 

-1.0493 

0.0147 

0.2235 

-0.003 

SF181 

04/07/03 

79.5 

356.2 

75 

700 

13 

12 

0 

0 

0 

17.4979 

0.1598 

-0.6601 

0,0119 

0.1675 

0.008 

SF182 

04/07/03 

79.5 

210.9 

100 

700 

10 

12 

0 

0 

0 

14.9671 

-0.0355 

-0.5625 

0.0014 

0.1805 

-0.0027 

SF183 

04/07/03 

79.5 

210.9 

75 

700 

13 

12 

0 

0' 0 

16.0427 

-0.0323 

-0.2541 

0.003 

0.1056 

-0.0015 

SF184 

04/07/03 

79.5 

501.5 

100 

700 

10 

12 

0 

0 

0 

2.9632 

-0.1106 

-0.8265 

0.0007 

0.2131 

-0.0102 

SF185 

04/07/03 

795 

501.5 

75 

700 

13 

12 

0 

0 

0 

0.0831 

-0.2192 

-0.3446 

-0.005 

0.1249 

-0.0096 







p 



Appendix G. AnalyzemodXls.m 

% AnalyzemodXls.m 
% Erik Oiler, 2003 

% Series Specific Treatment 

% Series 7: Inertial Results are not subtracted for Mass Matrix Tests. 
% Variables 

% ActualPhasefl ,f2,f3: The actual phase of the forces at the 
% oscillation frequency. 

% Ampl: The amplitude of the data point. 

% avg(l 3): The mean value of the force and motions for each 

% experiment. 

% B(6,6): The inverse sensitivity matrix for the load cell. 

% C(6): The values of capacitances in the electrical filters in 

% microfarads. 

% CF(6,6): The conversion factor matrix for the load cell. 

% CFtemp(6): The conversion factor vector for the load cell data. 
% Data(NuniPoints,15): The Datain matrix after interpolation for 
% constant time intervals. 


% Date: 

The date the experiment was performed. 

% DateNum: 

The numeric MATLAB code for the date and time the 

% experiment was 


% 

performed. 


% Datain; 

Holds all of the input data as read from the xls file. 

% 

Column 1: 

Y force 

% 

2 

Z Force 

% 

3 

X Force 

% 

4 

Pitch Moment 

% 

5 

Yaw Moment 

% 

6 

Roll Moment 

% 

7 

Reserved for wave data 

% 

8 

Reserved for wave data 

% 

9 

X position 

% 

10: Y position 

% 

1 

1: Z Position 

% 

12: Yaw Position 

% 

1. 

3: Pitch Position 

% 

14: Time 

% 

1 

5: Date 

% 

16: Time Interval between samples 

% DeltaT: 

The average time interval between samples in seconds. 

% Depth: 

The depth of the experiment in cm. 

% dropgap: 

The number of data samples dropped from the start of 


% the data in order to allow for acceleration. 

% DrivingFreq: The frequency of driving force oscillation in Hz. 

% This is omega converted into Hz. 

% Duration: The total length of time of the experiment. 

% existencecheck: Holds the results of the check for the input file. 

% fl: The ordered frequency of oscillation in Hz. 

% f2: The first harmonic of the oscillation frequency in Hz. 

% f3: The second harmonic of the oscillation frequency in Hz. 

% fname: The name of the file to be analyzed. 

% Fpitch: Calculated inertial force in pitch in Newton-meters. 
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% Fx: Calculated inertial force in the X direction in Newtons. 

% Fy: Calculated inertial force in the Y direction in Newtons. 

% Fyaw: Calculated inertial force in yaw in Newton-meters. 

% Fz: Calculated inertial force in the Y direction in Newtons. 

% Gain: The gain of the load cell amplifier. 

% infname: The name of the input file with the .xls extension added. 

% indexhigh: The index of the row of data representing the frequency 

% to begin filtering. 

% InertialForce: The Calculated Inertial Force in N or N-m. 

% Interval: The time duration of the fft data. 

% lyy: First moment of inertia abouthte Y-axis. 

% Izz: First moment of inertia about the Z-axis. 

% k: An index used in filtering. 

% LocalMax: The local maximum of data. Used in finding the peak 

% frequencies. 

% LocalMaxCounter: An index used to identify the lines of data for the 

% local maxima. 

% look: The size of the gap between frequency intervals. 

% manyrowsfid: The id number of the series row output file. 

% From the calling file. 

% mass: The mass of the model full of water in kg. 

% MotionPhi: The phase of the driving force motion at the frequency 

% of the local maxima. 

% nonzeroplanes: From the calling file. Indicates nonzero stem 
% planes or rudder. 

% NumLines: The number of lines in the input file containing force 
% and motion data. 

% NumLines2: Used to determine when the model stopped traveling in 
% the X direction. 

% NumPeriods: The number of periods that occur over the duration of 
% the experiment. 

% NumPoints: The number of points that will be used for the FFT. 

% omega: The frequency of oscillation in radians per second. 

% Period: The time in seconds of one cycle of motion. 

% PhaseIndexfl ,f2,f3: The line of Data containing the actual phase of 

% the motion at the oscillation frequency. 

% phi: The phase shift caused by dropping the first 1.2 seconds of 

% data. 

% Pitch Amp: Amplitude of Pitch oscillation in degrees. 

% PitchAngle: Steady pitch angle. Positive is nose up. 

% PitchDist: Ordered distance of travel in the Pitch direction in 

% degrees. 

% pitchfactor: The control factor for gantry motion in pitch. 

% PitchFreq: Frequency of Pitch oscillation in Hz. 

% PitchPhase: Ordered phase of Pitch oscillation in degrees. 

% PitchRadius: The distance from the center of gravity to the axis of 
% pitch in cm. 

% PitchSpd: Speed in the Pitch direction in degreees/sec. 

% PitchSupSpd: Superimposed speed in the Pitch direction in deg/sec. 

% PitchSupDir: Direction of superimposed motion. 

% R(6): The values of resistances in the electrical filters in 

% ohms. 

% Results{3,24): Holds all of the peaks of the FFT. 

% rowfname: The name of the row output file. 

% rowed: The file id number of the row output file. 

% Rudder: The angle of the rudder in degrees. 
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% series: The series of data to be analyzed. 

% Comes from the calling program. 

% 1: Horizontal Plane 

% 2: Vertical Plane 

% 3: Pure Sway 

% 4: Pure Heave 

% 5 : Pure Pitch 

% 6: Pure Yaw 

% 7: Mass Matrix 

% 8: Steady Force 

% 9: Check of Inertial Calculations 

% 10: Miscellaneous 

% 11: Oscillation Tests 

% SampleFreq: Calculated sample frequency based on average interval 
% between datapoints. 

% shiftlength: The distance between the center of gravity and midships 
% in m. 

% StemPlanes: The angle of the stem planes in degrees. 

% Submergence: The submergence of the model to the top of the hull in 
% mm. From the calling file. 

% Summary(3,24): Holds data from the FFT peak closest to the driving 
% frequency. 

% Time: Ordered length of time of the experiment in seconds. 

% TimeStart: The first recorded time of the experiment. 

% Vexc: The excitation voltage of the amplifier. 

% write: A switch for whether or not to write the individual row 
% output file. 

% XAmp: Amplitude of X oscillation in cm. 

% XDist: Ordered distance of travel in the X direction in cm. 

% xfactor: The control factor for gantry motion in the X direction. 

% XFreq: Frequency of X oscillation in Hz. 

% Xg: The distance from the model center of gravity to the 

% zero point of the load cell. 

% xhigh: The cutoff frequency for low pass filtering. 

% XPhase: Ordered phase of X oscillation in degrees. 

% XSpd: Speed in the X direction in cm/sec. 

% XSupSpd: Superimposed speed in the X direction in cm/sec. 

% XSupDir: Direction of superimposed motion. 

% YAmp: Amplitude of Y oscillation in cm. 

% Yaw Amp: Amplitude of Yaw oscillation in degrees. 

% Yaw Angle: The steady yaw angle in degrees. Positive is nose to 

% stbd. 

% YawDist: Ordered distance of travel in the Yaw direction in 


% degrees. 

% yawfactor: The control factor for gantry motion in yaw. 

% YawFreq: Frequency of Yaw oscillation in Hz. 

% YawPhase: Ordered phase of Yaw oscillation in degrees. 

% YawRadius: The distance from the center of gravity to the strut in 
% cm. 

% YawSpd: Speed in the Yaw direction in degreees/sec. 

% YawSupSpd: Superimposed speed in the Yaw direction in deg/sec. 

% YawSupDir: Direction of superimposed motion. 

% yfactor: The control factor for gantry motion in the Y direction. 

% YFreq: Frequency of Y oscillation in Hz. 

% YPhase: Ordered phase of Y oscillation in degrees. 

% YSpd: Ordered speed of travel in the Y direction in cm/sec. 
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% YSupSpd: 
% YSupDir: 
% ZAmp: 

% ZDist: 

% zfactor: 

% ZFreq: 

% ZPhase: 

% ZSpd: 

% ZSupSpd: 
% ZSupDir: 


Superimposed speed in the Y direction in cm/sec. 
Direction of superimposed motion. 

Amplitude of Z oscillation in mm. 

Ordered distance of travel in the Z direction in mm. 
The control factor for gantry motion in the Z direction. 
Frequency of Z oscillation in Hz. 

Ordered phase of Z oscillation in degrees. 

Ordered speed of travel in the Z direction in mm/sec. 

Superimposed speed in the Z direction in mm/sec. 
Direction of superimposed motion. 


% Initialize certain variables. 

NumPoints = 2048; 

Results = zeros(3,24); % Will hold all of the Peaks of the FFT. 

Summary = zeros(3,24); % Will hold the FFT peak closest to the driving 
% frequency. 

[Data] = 0; 

[Datain]=0; 


% Get the input data and find the length of the file, 
infname = strcat(fhame/.xls’); 

existencecheck = exist(infhame); % Determines if the input file exists, 
if existencecheck =0 

fprintf(’%s does not exist.\n',infname) 
return 
end 

Datain = xlsread(infhame); % Reads the input file. 

NumLines = size(Datain,l)-4; % Gets the number of data samples in the 
% input file. 


% Get the date from the input file and convert from EXCEL to 
% MATLAB format. 

Date=d atestr(Datain( l,15)-36525,2); 

% Build a date-time string to ensure the correct gains are applied. 
% Add 693960 to the date to get number of days from 0000. 

% Divide the number of minutes by 86400 to get fractional days. 
% MATLAB date serial numbers are in the form: 

% days since 0000.fraction of a day 

DateN um-Datain(1,15)+693960+Datain( 1,14)/86400; 


% Record input motion parameters. 
XAmp = Datain(NumLines+2,l); 
if isnan(XAmp)=0 

XFreq = Datain(NumLines+2,2); 
XPhase = Datain(NumLines+2,3); 
XSupSpd = Datain(NumLines+2,4); 
XSupDir = Datain(NumLines+2,5); 
else 

XAmp=0; 

XFreq-0; 

XPhase = 0; 

XSupSpd = 0; 

XSupDir = 0; 
end 


YAmp = Datain(NumLines+2,6); 




if isnan(YATnp)=0 

YFreq = Datain(NumLines+2,7); 
YPhase = Datain(NumLines+2,8); 
YSupSpd = Datam(NumLines+2,9); 
YSupDir = Datain(NumLines+2,l 0); 
else 

YAmp=0; 

YFreq = 0; 

YPhase = 0; 

YSupSpd = 0; 

YSupDir = 0; 
end 

ZAmp = Datain(NuniLines+2,l 1); 
if isnan(ZAmp)—0 

ZFreq = Datain(NumLines+2,12); 
ZPhase = Datam(NumLines+2,13); 
ZSupSpd = Datam(NuniLines+2,14); 
ZSupDir = Datain(NumLines+2,l 5); 
else 

ZAmp=0; 

ZFreq = 0; 

ZPhase = 0; 

ZSupSpd = 0; 

ZSupDir = 0; 
end 

YawAmp = Datain(NumLines+2,16); 
if isnan(YawAmp)=0 

YawFreq = Datain(NumLines+3,l); 
YawPhase = Datain(NuinLines+3,2); 
YawSupSpd = Datain(NumLines+3,3); 
YawSupDir = Datain(NumLines+3,4); 
else 

YawATnp=0; 

YawFreq = 0; 

YawPhase = 0; 

YawSupSpd = 0; 

YawSupDir = 0; 
end 

PitchAmp = Datain(NuniLines+3,5); 
if isnan(PitchAmp)—0 
PitchFreq = Datain(NumLines+3,6); 
PitchPhase = Datain(NumLines+3,7); 
PitchSupSpd = Datain(NumLines+3,8) 
PitchSupDir = Datain(NumLines+3,9); 
else 

PitchAmp^O; 

PitchFreq == 0; 

PitchPhase = 0; 

PitchSupSpd = 0; 

PitchSupDir = 0; 
end 


XSpd = DatainfNumLines+3,10); 




if isnan(XSpd)—1 
XSpd = 0; 

if ((isnan(XSupSpd)===0)&(XSupDir^=l)) 

XSpd = XSupSpd; 
end 
end 

XDist = Datain(NumLines+3,l 1); 
if isnan(XDist)=l 
XDist = 0; 
end 

YSpd = Datain(NuniLines+3,12); 

YDist = Datain(NumLines+3,13); 

ZSpd = Datain(NumLines+3,14); 

ZDist = Datain(NumLines+3,15); 

YawSpd = Datain(NumLines+3,16); 

YawDist = Datain(NumLines+4,l); 

PitchSpd = Datain(NumLines+4,2); 

PitchDist = Datain(NuinLines+4,3); 

Time = Datain(NumLines+4,5); 

% Determine actual sample frequency based upon the measured interval 
% between data points. Necessary because the control software did not 
% always sample at the ordered sample rate. 

SampleFreq = l/mean(Datain(l:NumLines,16)); 

% Determine the oscillation frequency. 
fl=0; 

0 = 0 ; 

B=0; 

if (XFreq ~ 0)&(isnan(XFreq)==0) 
fl = XFreq; 

0 = 2* XFreq; 

B = 3* XFreq; 
end 

if ((YFreq -= 0)&(isnan(YFreq)==0)) 
fl=YFreq; 

0 = 2* YFreq; 

B = 3* YFreq; 
end 

if (ZFreq ^ 0)&(isnan(ZFreq)=0) 
fl = ZFreq; 

0 = 2* ZFreq; 

B = 3* ZFreq; 
end 

if (PitchFreq — 0)&(isnan(PitchFreq)=0) 
fl =PitchFreq; 

0 = 2* PitchFreq; 

B = 3* PitchFreq; 




end 


# 



# 


t 










if (YawFreq — 0)&(isnan(YawFreq)=0) 
fl = YawFreq; 
f2 = 2* YawFreq; 

0 = 3* YawFreq; 
end 

omega=fl *2*pi; % The frequency of oscillation in radians per second. 

% Drop the first 1.2 seconds of data to account for acceleration. At 
% the sample frequency of 25 Hz, drop the first 30 points. Shift the 
% other points up in the array. 

dropgap = round(l .2*SampleFreq); 
for I = 1 iNumLines-dropgap; 

Datain(I,l :16)=Datain(I+dropgap,l :16); 
end 

NumLines = NumLines - dropgap; 

% Drop data recorded after the model stopped moving along the X-axis 

% for non-Inertial tests. 

if(XSpd-=0) 

NumLines2=NumLines; 
fori =NumLines-l:-l:10 

if (Datain(I,9)=Datain(I+l,9)) % Looks for constant X position. 

NumLines2 = I; 
end 
end 

NumLines = NumLines2; 
end 

% Make the time of good data be an integer number of wavelengths. 

Period =l/fl; 

Duration = Datain(NumLines, 14)-Datain( 1,14); 

NumPeriods = Duration/Period; 

NumLines = round(SampleFreq*Period*floor(NumPeriods)); 

% PREPROCESSING STAGE 
% Preprocess the input file for analysis. This includes: 

% 1) Accounting for the phase shift caused by dropping the first 1.2 
% seconds. 

% 2) Determining the mean force for each column and removing that mean 
% from the column data to get dynamic response, 

% 3) Converting force data from voltages to forces and moments, 

% 4) Converting position data to centimeters for x and y, millimeters 
% for z, and degrees for the angles, 

% 5) Converting the time past midnight to time of data run, 

% 6) Finding the Mean DeltaT, 

% 7) Interpolating the data and time to produce even time intervals 
% and the associated data. 

% 1,2) Account for the phase shift and determine the mean of the forces 
% columns. 

for J= 1:6 
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avg(J) = inean(Datain(l :NuinLines,J)); 
for I = liNumLines 

Datain(l,J) = Datain(l,J) - avg(J); 
end 
end 

% 3,4) Convert force data from voltages to forces and moments. 

% Convert position data to centimeters for x and y, millimeters for 
% z, and degrees for the angles, 

TimeStart = Datain(l,14); % Used for converting from time past midnight 
% to time of run. 

% The inverse sensitivity matrix (B) 

if ((datenum(Date)>=datenum(2002,07,17))&(datenum(Date)<datenum(2002,10,31))) 
[B] = [ 0.3901 0.0029 -0.0071 -0.0010 -0.0024 -0.0016; 

0.0023 0.3887 0.0016 0.0025-0.0039 0.0019; 

0.0139 0.0129 1.5006-0.0002-0.0215-0.0018; 

0.0000-0.0001 0.0016 0.0069 0.0000-0.0001; 

-0.0001 0.0000 0.0018 0.0000 0.0070 0.0000; 

-0.0003-0.0002 0.0000 0.0000 0.0000 0.0108]; 
else %valid on and after Oct 31, 2002 

[B] = [ 0.3856 O.OOlZO-0.0016-0.0008 0.0001 -0.0030; 

0.0023 0.3811 -0.0040 0.0012-0.0034 0.0031; 

0.0093 0.0024 1.5109 0.0068-0.0231 -0.0014; 

0.0001 0.0000 0.0013 0.0069 0.0000 0.0000; 

0.0000 0.0000 0.0007 0.0000 0.0069-0.0001; 

-0.0003 0.0001 0.0000 0.0000 0.0000 0.0106]; 

end 

% Establish Gains and Excitation Voltage 
if (DateNum>731434.61458) % August 6, 2002 1445, 

[Gain] = [4000;4000;1000;1000;1000;4000]; % y z x pitch yaw roll 
[Vexc] = [10; 10; 2.5; 10; 10; 10]; 
end 

% Calculate the Conversion Factors (CF) 

CFtemp = Gain.*Vexc*10^-6; 

CF = zeros(6,6); 

CF(l,l)=CFtemp(l); 

CF(2,2)= CFtemp(2); 

CF(3,3)= CFtemp(3); 

CF(4,4)= CFtemp(4); 

CF(5,5)= CFtemp(5); 

CF(6,6)= CFtemp(6); 

% Establish multipliers 

% yfactor has a negative multipler to convert from the gantry y positive 

% to the load cell y positive direction. 

if ((DateNum>731513.33333)&(DateNum<731551.33333)) 

% Between Oct 24,0800 and Dec 1,0800. 

xfactor = 3850; 

yfactor --2410; 

zfactor = 2114; 

pitchfactor = 972; 





yawfactor= 1818; 

elseifDateNum>731551.33333 %Dec 1,0800 
xfactor = 3850; 
yfactor =-2410; 
zfactor = 2114; 
pitchfactor= 155; 
yawfactor= 1818; 
else 

xfactor = 3850; 
yfactor =-2410; 
zfactor = 4921; 
pitchfactor = 1111; 
yawfactor= 1025; 
end 

% Convert voltages to forces and position data to the metric system, 
for 1 = 1 : NumLines 

Datain(I,l :6) = (CF^-1 *B*Datain(I,l :6yy; 

Datain(I,9) =Datain(l,9) /xfactor; 

Datain(I,10) = Datain(l,10) / yfactor; 

Datain(I,l 1) = Datain(l,l 1) / zfactor; 

Datain(l,12) = Datain(l,12) / yawfactor; 

Datain(I,13) = Datain(I,13) / pitchfactor; 

Datain(1,14) = Datain(l, 14) - TimeStart; 
end 

% Convert the position data to position of midships vice position of 
% the strut. 

distance = .1093; % The distance from the strut to midships. 

PitchRadius = .2576; % Length of pitch arm 
for 1 - 1: NumLines 

% Effect of yaw on y position. 

Datain(I,10) = Datain(I,10)-distance*sin(Datain(I,12)*pi/180); 

% Effect of pitch on x position. 

Datain(I,9) = Datain(l,9) + PitchRadius* sin(Datain(I,l l)*pi/180); 

% Effect of pitch on z position. 

Datain(I,l 1) =Datain(l,l l)-PitchRadius*(l-cos(Datain(l,l l)*pi/180)); 
end 

% Account for the transducer being forward of midships. 

% Reference: Marine Hydrodynamics by J.N. Newman, Sect 6.2 
% waveheight = Acos(kx-wt+arbitrary phase shift) (eq 7 page 240) 

% Let X = 0 at the midships location. 

% The wave is travelling in the positive x direction. 

% The transducer is at a negative x direction from midships. 

% Assume the wave is travelling in the positive x-direction. 

% Use midships as the zero reference location. 

% Shift the origin of the coordinate system from the origin of the load 
% cell to vessel amidships, 
shiftlength - 0.0522; 

wavenum = omega^2/9.81; 
for I =l:NumLines 

Datain(l,4) = Datain(I,4) - Datain(I,2) * shiftlength; % pitch 
Datain(I,5) = Datain(I,5) + Datain(I,l) * shiftlength; % yaw 




% 6) Determine the Delta! from the input file. 

Delta! == 1/SampleFreq; 

% 7) Interpolate data points to produce even time intervals and 
% associated data. By this process the input datain matrix will be 
% interpolated into the data matrix. 

% Start by building even time intervals and the date column, 
for I = 1 iNumPoints; 

% Interpolate as if the sample freq was 25 hz. 

Data(I,14) = (l-l)* .04; 

Data(1,15) = Datain( 1,15); 
end 

% Interpolate each column in turn, 
for J= 1:13 

Data(l,J) = Datain(l,J); 
for I = 2:NumLines; 
for K = 1 rNumLines; 

if(Datain(K,14)>Data(I,14)) 

break 

end 

end 

KS =K - 1; 

if (Data(1,14)<=Datain(NumLines, 14)) 

Data(I,J) - Datain(KS,J) + (Data(I,14)-Datain(KS,14))*... 
(Datain(K,J)-Datain(KS,J))/(Datain(K,l 4)-Datain(KS,l 4)); 

else 

Data(!J) = 0.0; 
end 
end 
end 


% Identify the driving force, 
if YFreq 0 & isnan(YFreq)=0 
DrivingForce = 10; 
omega = YFreq*2*pi; 
elseif ZFreq— 0 & isnan(ZFreq)==0 
DrivingForce = 11; 
omega = ZFreq*2*pi; 
elseif YawFreq— 0 & isnan(YawFreq)==0 
DrivingForce = 12; 
omega = YawFreq*2*pi; 
elseif PitchFreq~0 & isnan(PitchFreq)==0 
DrivingForce = 13 ; 
omega = PitchFreq*2*pi; 
elseif XFreq ~ 0 & isnan(XFreq)==0 
DrivingForce = 9; 
omega = XFreq*2*pi; 
else 

DrivingForce=9; 

omega=0; 




% PROCESSING STAGE 
for J = 10:13 

Data(l :NumLines,J)=Data(l :NumLines,J)-mean(Data(l :NumLines,J)); 
end 


% Round all data up to 2048 points. Pad the excess with zero’s. 
Data(NumLines+l :NumPoints,l :13) = 0.0; 

% Convert the time column to a frequency column. 

Interval = .04 * NumPoints; 

for I = 1: NumPoints 
Data(I,l 4)=(I-1 yinterval; 
end 

% For each column, perform Fast Fourier Transformation, low pass 
% filtering, and identify the frequencies and amplitudes of the 
% forces. 

% Perform the Fast Fourier Transformation. 
forJ=l: 13 

Data(:,J) = fft(Data(:,J)); 

% Transform the FFT Coefficient to account for padding, 
for I = 1 -.NumPoints 

Data(I,J) = Data(l,J)* NumPoints/NumLines; 
end 
end 


% Change the phase of all data such that the phase is relative to 
% the driving motion. 


for I = 1 :NumLines 
for 1:13 

Data(I,J) = Data(l,J)*exp{-i*angle(Data(I,DrivingForce))); 
end 
end 

% Account for ELECTRICAL FILTERS. 

[R] = [145500 144100 147400 149200 191600 146600]; 

[C] = [.073 .114 .081 .099 .089 .075 ]*10^-6; 

for J= 1:6 

for I = 1 iNumLines 

Data(I,J) = Data(I,J)*exp(-i*atan(omega*C(J)*R(J))) ; 
end 
end 


% Determine the actual phase of the motion at the oscillation frequency 
% Find the index of the frequency of oscillation. 
look= .5*1/Interval; 

Phaselndexfl == find(((fl dook)<=Data(:,14))&(Data(:,14)<(fl+look))); 
PhaseIndexfZ = find(((f2-look)<=Data(:, 14))&(Data(:, 14)<(f2+look))); 



# 


PhaseIndexO = find(((f3-look)<=Data(:,l 4))&(Data{:,14)<(B+look))); 
for J=9:13 

ActualFreq( J)=Data(PhaseIndexfl ,14); 

H ActualPhasefl (J)= angle(Data(PhaseIndexfl ,J))* 180/pi; 

ActualPhasef2(J)" angle(Data(Phaselndexf2,J))*l 80/pi; 
ActualPhasef3(J)= angle(Data(PhaselndexO,J))* 180/pi; 
ActualMotion(J)=sqrt(2.*2.*((abs(Data(Phaselndexfl ,J))/NumPoints)'^2)); 
end 

Jk omega = ActualFreq(DrivingForce)*2*pi; 









% Find forces at the frequencies of oscillation, 
for J == 1:6 

Summaryfl(3*J-2)=Data(PhaseIndexfl,14); % freq 

Summaryfl(3*J-l)=sqrt(2.*2.*((abs(Data(PhaseIndexfl ,J))/NuniPoints)^2)); % ampl 
Summaryfl (3*J)=angle(Data(PhaseIndexfl ,J))* 180/pi-ActualPhasefl (DrivingForce); %phase 

Summaryf2(3*J-2)=Data(PhaseIndexf2,14); % freq 

Summaryf2(3*J-l)=sqrt(2.*2.*((abs(Data(PhaseIndexf2,J))/NumPoints)'^2)); % ampl 
Summaryf2(3*J)=angle(Data(PhaseIndexf2,J))* 180/pi-ActualPhasef2(DrivingForce); %phase 

SummaryB(3 * J-2)=Data(Phaselndexf3,14);% freq 

Summaryf3(3*J-l)=sqrt(2.*2.*({abs(Data(PhaseIndexf3,J))/NumPoints)^2));% ampl 
Summaryf3(3*J)=angle(Data(Phaselndexf3,J))* 180/pi-ActualPhasef3(DrivingForce); %phase 
end 

% Calculate the inertial forces in mks units for all series except 

% mass matrix evaluation (7). 

mass=3.67; 

Xg=0.03; 

Gyradius = 0.0764; 

Ya^vRadius = 0.07; 

if series ~ 7 % Calculates inertial force for all except mass matrix 
% (inertial) tests, 
lyy = 0.092509; 

Izz = 0.1006; 

XAmp - ActualMotion(9); 

YAmp = ActualMotion(lO); 

ZAmp = ActualMotion(l 1); 

PitchAmp = ActualMotion(13); 

YawAmp = ActualMotion(12); 

Fx = -omega^2*mass*(XAmp/100 + PitchRadius* (Pitch Amp*pi/180)); 

Fy = -omega^2'*'(mass*(YAmp/100 + (Yavv^Amp*pi/180)*YawRadius) + mass *Xg* (Yaw Amp *pi/l 80)); 
Fz - ♦‘Omega^2*(mass*(ZAmp/l 000 + (PitchAmp*pi/l 80) * (YawRadius))- 
mass*Xg*(PitchAmp*pi/l 80)); 

Fpitch = -omega^2*(-mass*Xg*(ZAmp/1000 + (PitchAmp*pi/180) * (YawRadius-Xg)) + lyy * 
(PitchAmp*pi/l 80)); 

Fyaw = ‘Omega^2 * (mass*Xg*(YAmp/100 - (YawAmp*pi/180)*YawRadius)4-Izz*(YawAmp*pi/l 80)); 
else 

lyy = 0; 
lzz=0; 

Fx=0; 

Fy=0; 

Fz=0; 

Fpitch = 0; 
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% Calculate the hydrodynamic response at the oscillation frequency and 
% its harmonics. 

DrivingFreq = omega/(2*pi); 


forJ = l:6 
switch J 
case 1 

InertialForce = Fy; 
case 2 

InertialForce = Fz; 
case 3 

InertialForce = Fx; 
case 4 

InertialForce = Fpitch; 
case 5 

InertialForce = Fyaw; 
case 6 

InertialForce == 0; 
end 


Results(l,4*J’3) = Summaryfl(3*J-2); % oscillation freq 
Results(l,4*J-2) = Summaryfl(3*J-l); % total force 

Results(l,4*J-l) = imag(Summaryfl(3*J-l))+real(Summaryfl(3*J-i)) + InertialForce; % force without 
inertial component 

Results(l,4*J) = Summaryfl(3*J); % phase 


Results(2,4*J-3) = Summaryf2(3*J-2); % first harmonic of oscillation freq 
Results(2,4*J-2) = Summaryf2(3*J-l); % total force 
Results(2,4*J-l) = Summaryf2(3*J-l);% 

Results(2,4*J) = Summaryf2{3*J); % phase 


Results(3,4*J-3) = Summary£3(3*J-2); % second harmonic of oscillation freq 
Results(3,4*J-2) = Summaryf3(3% total force 
Results(3,4*J-l) = Summaryf3(3*J“l); % force without inertial component 
Results(3,4*J) = Summaryf3(3*J); % phase 

end 

% Condition the Results array. 

% Ensure the phase angles are +-180 degrees. 

% Ensure the amplitudes are all positive, 
for I =1:3 
for J = 1:6 

if Results(I,4*J-l) <0 

Results(I,4*J-l) = abs(Results(I,4*J-l)); 

Results(I,4*J) = Results(I,4*J) + 180; 
end 

ifResults(I,4*J)>180 
Results(I,4*J) = Results(I,4*J) - 360; 
elseif Results(l,4*J) <-180 

Results(I,4*J) = Results(I,4*J) + 360; 




% Print to a common row output file. 
fprintf(manyrowsfid,*%s \t %s\t %6.2f,fname,Date,Depth); 

fyrintf(manyrowsfid,V/o6.4f \t%5.1f \t %5.2f\t%5. If,Submergence,XSpd,XDist,Duration); 
fprintf(manyrowsfid,’ \t %2.0f \t%2.0ft %3.1f.t %3.1f ,Pitch Angle,Yaw Angle,StemPlanes,Rudder); 
fprintf(manyrowsfid,’ \t %8.5f \t %7,4f \t %4.1f \t %4.1f \t%8.5f \t %7.4f \t %7.4f \t 
%4.1 f ,ActualFreq(9),ActualMotion(9),XPhase, ActualPhasefl (9),Results( 1,9:12)); 
fprintf(manyrowsfid,’ \t %8.5f \t %7.4f\t %7.4f \t %4.1f,Results(2, 9:12)); 
fprintf(manyrowsfid,’ \t %8.5f \t %7.4f \t %7.4f \t %4.1f,Results(3,9:12)); 
fprintf(manyrowsfid,' \t %8.5f\t %7.4f\t %4.1f\t %4.1f\t %8.5f \t %7.4f \t %7.4f\t 
%4.1 f ,ActualFreq( 10), ActualMotion( 10),YPhase, ActualPhasefl (10),Results( 1,1:4)); 
fprintf(manyrowsfid,' \t %8.5f \t %7.4f\t %7.4f\t %4.1f,Results(2,1:4)); 
fprintf(manyrowsfid,’ \t %8.5f \t %7.4f \t %7.4f \t %4.1f,Results(3,1:4)); 
fprintf(manyrowsfid,’ \t %8.5f \t %7.4f \t %4,lf \t %4.1f \t %8.5f \t %7.4f \t %7.4f \t 
%4.1f,ActualFreq(l l),ActualMotion(l l),ZPhase, ActualPhasefl (1 l),Results(l,5:8)); 
^rintf(manyrowsfid,’ \t %8.5f\t %7.4f \t %7.4f \t %4.1f,Results(2, 5:8)); 
fprintf(manyrowsfid,’ \t %8.5f\t %7.4f\t %7.4f \t %4.1f,Results(3, 5:8)); 
fprintf(manyrowsfid,’ \t %8.5f \t %7.4f \t %7.4f \t %4.1f,Results(l,21:24)); 
fprintf(manyrowsfid,’ \t %8.5f\t %7.4f \t %7.4f\t %4.1f,Results(2, 21:24)); 
fprintf(manyrowsfid,’ \t %8.5f\t %7.4f \t %7.4f\t %4.1f,Results(3,21:24)); 
fprintf(manyrowsfid,’ \t %8.5f \t %7.4f \t %4.1f \t %4.1f \t %8.5f \t %7.4f \t %7.4f \t 
%4.1 f, ActualFreq( 13),ActualMotion( 13) ,PitchPhase, ActualPhasefl (13),Results( 1,13:16)); 
fprintf(manyrowsfid,’ \t %8.5f \t %7.4f \t %7.4f \t %4.1f,Results(2,13:16)); 
fi)rintf(manyrowsfid,’ \t %8.5f \t %7.4f \t %7.4f \t %4.1f,Results(3,13:16)); 
fi>rintf(manyrowsfid,’ \t %8.5f \t %7.4f\t %4.1f\t %4.1f \t %8.5f \t %7.4f \t %7.4f\t 
%4.1 f, ActualFreq( 12), ActualMotion( 12), Y awPhase, ActualPhasefl (12),Results( 1,17:20)); 
fprintf(manyrowsfid,’ \t %8.5f\t %7.4f \t %7.4f\t %4.1f,Results(2, 17:20)); 
fprintf(manyrowsfid,' \t %8.5f \t %7.4f \t %7.4f \t %4.1f \n',Results(3, 17:20)); 




Appendix H. Results of Inertial Force Calculation Checks 












































Appendix L Selected Portion of the Output File from AnalyzemodXls.m for the Analysis of 
the Oscillation Test Series 

The first two rows have been added by the author for clarity._ 
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Appendix J. CoeffSolver.m 

%CoeffSolver.in 


% Erik Oiler 
% Spring 2003 

% Solves for coefficients using linear regression. 

% Added mass and damping coefficients are of the form 
% Coeff = aO + al *L/subm + a3*Fr 
% Restoring forces equations are of the form 
% Coefficient=aO + al*L/subm + a2*L/subm^2 +a3*Fr 
clear all; 
close all; 


% Solve for Restoring Force Coefficients 

% Solve for Yaw Induced Restoring Force (Yuv) 

YuvMatrix = [-5.368865729 

-4.853454974 

-2.329293369 

-3.514754719 

-2.848702047 

-1.614922323]; 


SteadyY awMatrix=[ 1 
1 1.720624859 

1 1.260063003 

1 2.711810649 

1 1.720624859 

1 1.260063003 


2.711810649 

2.960549904 

1.587758772 

7.353916996 

2.960549904 

1.587758772 


7-353916996 0.289533426 

0.289533426 

0.289533426 

0.386044568 

0.386044568 

0.386044568]; 


YuvCoeffMatrix=SteadyYawMatrixVYuvMatrix; 

PredY uvMatrix=Steady Y awMatrix* Y uvCoeffMatrix; 
for i = 1 :size(PredYuvMatrix,l) 

YuvDiff(i) = abs((PredYuvMatrix(i)-YuvMatrix(i)))A^uvMatrix(i); 
end 

Y uvrms=norm(Y uvDiff)/sqrt(size( Y u vMatrix, 1)); 


% Solve for Pitch Induced Restoring Force (Zuw), C53 

ZuwMatrix = [1.023955274 

2.269136274 

1.141877357 

1.220249361 

2.612356757 

1.801033337]; 


SteadyPitchMatrix=[l 

2.711810649 

7.353916996 0.289533426 

1 

1.720624859 

2.960549904 

0.289533426 

1 

1.260063003 

1.587758772 

0.289533426 

1 

2,714285714 

7.367346939 

0.386044568 

1 

1.718592965 

2.953561779 

0.386044568 

1 

1.259668508 

1.586764751 

0.386044568]; 


ZuwCoeffMatrix=SteadyPitchMatrix\Zu wM atrix; 






PredZuwMatrix=SteadyPitchMatrix*ZuwCoeffMatrix; 
for i = 1 ;size(PredZuwMatrix,l) 

ZuwDiff(i) = abs((PredZuwMatrix(i)-ZuwMatrix(i)))/ZuwMatrix(i); 
^ end 

Zuwnns=norm(ZuwDiff)/sqrt(size(ZuwMatrix, 1)); 


% Solve for Yaw Induced Restoring Moment (Nuv) 


NuvMatrix = [-1.020609587 

-0.996711818 

-0.682640735 

-0.800283816 

-0.737734321 

-0.586442502]; 


SteadyY awMomMatrix=[ 1 
1 1.720624859 

1 1.260063003 

1 2.711810649 

1 1.720624859 

1 1.260063003 


2.711810649 7.353916996 


2.960549904 

1.587758772 

7.353916996 

2.960549904 

1.587758772 


0.289533426 

0.289533426 

0.386044568 

0.386044568 

0.386044568]; 


0.289533426 


^ NuvCoeffMatrix=SteadyYawMomMatrix\NuvMatrix; 

PredNuvMatrix=SteadyYawMomMatrix*NuvCoeffMatrix; 
for i = 1 :size(PredNuvMatrix,l) 

NuvDiff(i) = abs((PredNuvMatrix(i)-NuvMatrix(i)))/NuvMatrix(i); 
end 

Nuvnns=norm(NuvDifi)/sqrt(size(NuvMatrix, 1)); 


% Solve for Pitch Induced Restoring Moment (Muw) 
MuwMatrix = [-0.626004966 


-0.594465093 

-0.71149292 

-0.401198364 

-0.619156216 

-0.599156406]; 


SteadyPitchMomMatrix=[ 1 


2.711810649 7.353916996 


1.720624859 

1.260063003 

2.714285714 

1.718592965 

1.259668508 


2.960549904 

1.587758772 

7.367346939 

2.953561779 

1.586764751 


0.289533426 

0.289533426 

0.386044568 

0.386044568 

0.386044568]; 


0.289533426 


MuwCoeffMatrix=SteadyPitchMomMatrix\MuwMatrix; 

PredMuwMatrix=SteadyPitchMomMatrix*MuwCoeflMatrix; 

® for i = l:size(PredMuwMatrix,l) 

MuwDiff(i) = abs((PredMuwMatrix(i)-MuwMatrix(i)))/MuwMatrix(i); 
end 

Muwrms=norm(MuwDiff)/sqrt(size(MuwMatrix,l)); 


% Solve for Sway induced added masses 


% Solve for Yvdot 
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YvdotMatrix=[-0.017470 

-0.018969 

-0.016666 

-0.019724 

-0.017526]; 


SwayMatrix=[ 1 1.277450258 0.127678512 

1 2.749405234 0.383418954 

1 1.744402516 0.255740442 

1 1.277450258 0.383418954 

1 2.749405234 0.127678512]; 

YvdotCoeffMatrix = SwayMatrixWvdotMatrix; 
PredY vdotMatrix=SwayMatrix* Y vdotCoeffMatrix; 


sumofsquares=0; 


for i = l:size(PredYvdotMatrix,l) 

YvdotDiff(i) = (PredYvdotMalTix(i)-YvdotMatrix(i))A'vdotMatrix(i); 
end 

Yvdotrms=norm(YvdotDiflf)/sqrt(size(YvdotMatrix,l)); 

% Solve for Yv 

YvMatrix=[-0.058883 

-0.027836 

-0.062066 

-0.027142 

-0.058392]; 

YvCoeffMatrix = SwayMatrixVYvMatrix; 

PredY vMatrix=SwayMatrix* Y vCoeffMatrix; 

sumofsquares=0; 

for i = 1: size(Pred Y vMatrix, 1) 

YvDiff(i) - abs((PredYvMatrix(i)-YvMatrix(i)))A'vMatrix(i); 
end 

Y vnns=nonn(YvDifE)/sqrt(size(Y vMatrix, 1)); 


% Solve for Nvdot 

NvdotMatrix=[-0.000809 

-0.000692 

-0.000677 

-0.000862 

-0.000833]; 

NvdotCoeffMatrix = SwayMatrix\NvdotMatrix; 
PredNvdotMatrix=SwayMatrix*NvdotCoeffMatrix; 

for i = 1 :size(PredNvdotMatrix, 1) 

NvdotDiff(i)=abs(PredNvdotMatrix(i)-NvdotMatrix(i))/NvdotMatrix(i); 

end 

Nvdotrms=norm(NvdotDiff)/sqrt(size(NvdotMatrix,l)); 

% Solve for Nv 

NvMatrix=[-0.004323 

-0.008228 





-0.005161 

-0.008666 

-0.004136]; 

NvCoeffMatrix = SwayMatrix\NvMatrix; 
PredNvMatrix=SwayMatrix*NvCoefflVlatrix; 

for i = 1 :size(PredNvMatrix,l) 

NvDiff(i) = (PredNvMalTix(i)-NvMatrix(i))/NvMatrix(i); 

end 

Nvrms=norm(NvDiff)/sqrt(size(NvMatrix,l)); 

% Solve for coefficients in heave 
% Solve for Zwdot 

ZwdotMatrix=[-0.019162 
-0.022988 
-0.016744 
-0.020764 
-0.021686]; 


HeaveMatrix=[l 

1.421 

1 

2.547 

0.383 

1 

1.744 

0.256 

1 

1.421 

0.383 

1 

2.547 

0.128]; 


ZwdotCoeffMatrix = HeaveMatrixVZwdotMatrix; 
PredZwdotMatrix=HeaveMatrix*ZwdotCoeffMatrix; 

for i = 1 :size(PredZwdotMatrix,l) 

ZwdotDiff(i)=abs((PredZwdotMatrix(i)-ZwdotMatrix(i)))/ZwdotMatrix(i); 

end 

Zwdotrms=nonn(ZwdotDiff)/sqrt(size(ZwdotMatrix, 1)); 

% Solve for Zw 

ZwMatrix=[-0.062616 

-0.029098 

-0.065158 

-0.024766 

-0.063094]; 

ZwCoeffMatrix = HeaveMatrixXZwMatrix; 
PredZwMatrix=HeaveMatrix*ZwCoeffMatrix; 

for i = 1 :size(PredZwMatrix,l) 

ZwDiff(i) = abs((PredZwMatrix(i)-ZwMatrix(i)))/ZwMatrix(i); 
end 

Zv^rms=norm(Zv/Diff)/sqrt(size(ZwMatrix, 1)); 


% Solve for Mwdot 

MwdotMatrix=[0.002835 

0.002002 

0.002847 




0.000781 

0.003218]; 

MwdotCoeffMatrix = HeaveMatrix\MwdotMatrix; 
PredMwdotMatrix=HeaveMatrix*MwdotCoeffMatrix; 

for i = 1 :size(PredMwdotMatrix,l) 

MwdotDiff(i) = abs((PredMwdotMatrix(i)"MwdotMatrix(i)))/MwdotMatrix(i); 
end 

Mwdotnns=norm(MwdotDiff)/sqrt(size(MwdotMatrix,l)); 

% Solve for Mw 

MwMatrix=[0.018451 

0.008741 

0.009078 

0.010553 

0.012233]; 

MwCoeffMatrix = HeaveMatrix\MwMatrix; 
PredMwMatrix=HeaveMatrix*MwCoeffMatrix; 

for i = 1 :size(PredMwMatrix,l) 

MwDiff(i) = abs((PredMwMatrix(i)-MwMatrix(i)))/MwMatTix(i); 
end 

Mwnns=norm(MAvDiff)/sqrt(size(MwMatrix,l)); 


% Solve for coefficients in yaw 
% Solve for Yrdot 


YrdotMatrix=[0.001338 

0.002209 

0.003194 

0.002134 

0.001086]; 


YawMatrix-[l 
1 1.277 

1 1.744 

1 2.749 

1 1.277 


2.749 0.128 

0,128 

0.256 

0.383 

0.383]; 


YrdotCoeffMatrix = YawMatrix\YrdotMatrix; 

Pred YrdotMatrix=YawMatrix* YrdotCoeffMatrix; 

for i ^ l:size(PredYrdotMatrix,l) 

Y rdotDiff{i)=abs((PredYrdotMatrix(i)-YrdotMatrix(i)))/Y rdotMatrix(i); 
end 

Y rdotrms=norm(Y rdotDiff)/sqrt(size(YrdotMatrix, 1)); 


% Solve for Yr 

YrMatrix-[0.021030 

0.018478 

0.015644 

0.014240 





0.016104]; 


YrCoeffMatrix = YawMatrixWrMatrix; 
PredYrMatTix=Y awMatrix* Y rCoeffMatrix; 


for i = 1 :size(PredYrMatrix,l) 

YrDiff(i) = abs((PredYrMatrix(i)-YrMatrix(i)))/Y rMatrix(i); 
end 

YiTms=norm(YrDifO/sqrt(size(Y rMatrix, 1)); 

% Solve for Nrdot 

NrdotMatrix=["0.000871 

-0.000844 

-0.00113666 

-0.001051 

-0.000960]; 

NrdotCoeffMatrix = YawMatrix\NrdotMatrix; 
PredNrdotMatrix=YawMatrix*NrdotCoeffMatrix; 

for i = 1 :size(PredNrdotMatrix,l) 

NrdotDiff(i)=abs((PredNrdotMatrix(i)-NrdotMatrix(i)))/NrdotMatrix(i); 

end 

Nrdotrms=norm(NrdotDiff)/sqrt(size(NrdotMatrix,l)); 

% Solve for Nr 

NrMatrix=[-0.008514 

-0.005804 

-0.003476 

-0.004136 

-0.003916]; 

NrCoeffMatrix = Y awMatrixVNrMatrix; 
PredNrMatrix=YawMatrix*NrCoeffMatrix; 

for i = 1 :size(PredNrMatrix,l) 

NrDiff(i) = abs((PredNrMalTix(i)-NrMatrix(i)))/NrMatrix(i); 
end 

Nrrms=norm(NrDiff)/sqrt(size(NrMatrix,l)); 

% Solve for coefficients in Pitch 
% Solve for Zqdot 

ZqdotMatrix=[-0.001497196 

-0.000740791 

0.0036744]; 

PitchMatrix=[l 2.749405234 0.127678512 

1 1.277450258 0.127678512 

1 2.749405234 0.383418954]; 

ZqdotCoeffMatrix = PitchMatrix\ZqdotMatrix; 
PredZqdotMatrix=PitchMatrix*ZqdotCoeffMatrix; 
for i = 1 :size(PredZqdotMatrix,l) 

ZqdotDiff(i)=abs((PredZqdotMatrix(i)-ZqdotMatrix(i)))/ZqdotMatrix(i); 




end 

Zqdotrms=norm(ZqdotDiff)/sqrt(size(ZqdotMatrix,l)); 


^ % Solve for Zq 

ZqMatrix=[-0.043506132 

-0.048655759 

-0.009306293]; 



# 



ZqCoeffMatrix = PitchMatrix\ZqMatrix; 
PredZqMatrix=PitchMatrix*ZqCoeffMatrix; 

fori = l:size(PredZqMatrix,l) 

ZqDiff(i) = abs((PredZqMatrix(i)-ZqMatrix(i)))/ZqMatrix(i); 
end 

Zqrms=norm(ZqDiff)/sqrt(size(ZqMatrix, 1)); 

% Solve for Mqdot 
MqdotMatrix=[-0.000825618 
-0.00098228 
-0.003345578]; 

MqdotCoeffMatrix - PitchMatrix\MqdotMatrix; 
PredMqdotMatrix=PitchMatrix*MqdotCoeffMatrix; 

for i = 1 :size(PredMqdotMatrix, 1) 

MqdotDiff(i)=abs((PredMqdotMatrix(i)-MqdotMatrix(i)))/MqdotMatrix(i); 

end 

Mqdotrms=nonn(MqdotDiff)/sqrt(size(MqdotMatrix,l)); 


% Solve for Mq 

MqMatTix~[-0.008513669 

-0.007481107 

-0.000642713]; 

MqCoeffMatrix = PitchMatrixXMqMatrix; 

PredMqMatrix=PitchMatrix*MqCoeffMatrix; 

for i = 1 :size(PredMqMatrix,l) 

MqDiff(i) = abs((PredMqMatrix(i)-MqMatrix(i)))/MqMatTix(i); 
end 

Mqrms=nonn(MqDiff)/sqrt(size(MqMatrix, 1)); 

% Build Output File 
v^aming off 
delete('CoeffSolver.txf) 
w^aming on 

fileid = fopen(’CoeffSolver.txtVa'); 

fprintf(fileid;%s\t%s\t%s\t%s\t%s\t%s\t\n’;Coeff;r;L/Subm\’(L/Subm)^2’,*Fr’;rmsO; 
fi)rintf(fileid;%s\t%8.6f\t%8.6At%8.6f\t%8.6f\t%g\nVYuv’,YuvCoeffMatrix,Yuvrnis); 
fprintf(fileid/%s \t %8.6f \t %8.6f\t %8.6At %8.6f\t%g \nVZuw*,ZuwCoeffMatrix,Zuwrms); 
fprintf(fdeid;%s \t %8.6f \t %8.6fvt %8.6f \t %8.6f \t%g \n’;Nuv’,NuvCoeffMatrix,Nuvrms); 
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fprintf(fi!eid,'%s \t %8.6f \t %8.6f\t %8.6f \t%8.6f \t %g \nVMuw',MuwCoeffMatrix,Muwrms); 
fprintf(fileid;%s \t%s \t%s \t%s \t%s \t\n’,’Coeff;i ',’L/Subm’,’Fr',’nns’); 
fprintf(fileid,’%s \t %8.6f\t %8.6f\t %8.6f\t %g \n’,’Yrdot’,YrdotCoef{Matrix,Yrdotrms); 
fprintf(fileid,'%s \t %8.6f\t %8.6f\t %8.6At %g \n',Yr’,YrCoeffMatrix,Yrrms); 
fprintf(fileid,'%s \t %8.6f\t %8.6f\t %8.6At %g \n’,Yvdot’,YvdotCoeffMatrix,Yvdotnns); 
fprintf(fileid,'%s \t %8.6f\t %8.6f\t %8.6At %g \nVYv’,YvCoeffMatrix,Yvnns); 
fprintf(fileid,'%s \t %8.6f \t %8.6f\t %8.6f\t %g \n',’Zqdot’,ZqdotCoeffMalTix,Zqdotnns); 
fprintf(fileid,’%s \t %8.6f\t %8.6At %8.6i\t %g \nVZq’,ZqCoeffMatrix, Zqrms); 
Q)rintf(fileid,'%s \t %8.6f\t %8.6f\t %8.6f\t %g \n',’Zwdot’,ZwdotCoeffMatrix,Zwdotnns); 
fprintf(fileid,'%s \t %8.6f\t %8.6f\t %8.6f\t %g \n’,’Zw’,ZwCoeffMatrix, Zwrms); 
fprintf(fileid,'%s \t %8.6f\t %8.6f\t %8.6f\t %g \n’,'Mqdot’,MqdotCoeffMatrix,Mqdotnns); 
fprintf(fileid,'%s \t %8.6f\t %8.6f\t %8.6f\t %g \n’,’Mq’,MqCoeffMatrix, Mqrms); 
fprintf(fileid,’%s \t %8.6f\t %8.6f\t %8.6f\t %g \n’,’Mwdot’,MwdotCoeffMatrix,Mwdotnns); 
fprintf(fileid,'%s \t %8.6f \t %8.6f\t %8.6At %g \n’,’Mw’,MwCoefiOVIatrix, Mwrms); 
fprintf(fileid,'%s \t %8.6f\t %8.6f\t %8.6At %g \n','Nrdot',NrdotCoeffMatrix, Nrdotrms); 
l^rintf(fileid,’%s \t %8.6f\t %8.6f\t %8.6f\t %g \n’;Nr’,NrCoeffMatrix, Nrrms); 
fprintf(fileid,'%s \t %8.6f \t %8.6At %8.6f\t %g \nVNvdot’,NvdotCoeffMatrix, Nvdotrms); 
fprintf(fileid,'%s \t %8.6f\t %8.6f\t %8.6At %g \n',’Nv’,NvCoeffMatrix, Nvrms); 
status = fclose(fileid); 




Appendix K. Output of CoeffSolver.m 

This file shows the output of the file CoeffSolver.m. To derive the equation for 
coefficient A', find the sum of the products of the elements in the row labeled A and the 
value of the header for the column of the element, excepting the last column. In 
mathematical form, the equation looks is 

A' = y'{A-Header.) (44) 

y=i 

where n =4 for restoring force coefficients and n=3 for all other coefficients. The root 
mean square of the difference between the predicted and measured values are shown in 
the final column. 


Coeff 

1 

L/Subm 

(L/Subm)''2 

Fr 

rms 

Yuv 

3.035943 

-11.2325 

2.39970 

15.79519 

0.132812 

Zuw 

-7.723764 

9.156976 

-2.36537 

4.182974 

0.0752513 

Nuv 

0.020017 

-1.45297 

0.31796 

1.987686 

0.050452 

Muw 

-1.168727 

0.127681 

-0.00757 

1.079278 

0.0951314 

Coeff 

1 

L/Subm 

Fr 

rms 


Yrdot 

0.002324 

-8.7E-05 

-0.00063 

0.462378 


Yr 

0.020906 

0.000403 

-0.01798 

0.072234 


Yvdot 

-0.016345 

0.000061 

-0.00722 

0.043182 


Yv 

-0.081458 

0.001774 

0.12175 

0.148568 


Zqdot 

-0.002666 

-0.00051 

0.02022 

2.00E-16 


Zq 

-0.070199 

0.003498 

0.13373 

8.66E-16 


Zwdot 

-0.013633 

-0.00268 

-0.00567 

0.088254 


Zw 

-0.086685 

0.00091 

0.14075 

0.170412 


Mqdot 

0.00014 

0.000106 

-0.00985 

1.66E-16 


Mq 

-0.010515 

-0.0007 

0.03078 

1.18E-15 


Mwdot 

0.002825 

0.000594 

-0.00641 

0.271656 


Mw 

0.023224 

-0.00294 

-0.02236 

0.190629 


Nrdot 

-0.000786 

-1.9E-05 

-0.00058 

0.080731 


Nr 

-0.006013 

-0.00117 

0.01230 

0.232831 


Nvdot 

-0.00089 

0.000037 

0.00017 

0.092017 


Nv 

-0.00207 

0.000093 

-0.01649 

0.093355 










Appendix L. Model Scale Experiments Performed at MIT to 
Determine the Restoring Forces and Moments due to Body 
Angle 


Test # 

Water 

Depth 

Subm. 

Yaw 

Pitch 

Velocity 


m 

m 

Degrees 

Degrees 

m/sec 

SF98 

0.795 

0.3962 

0 

0 

1 

SF99 

0.795 

0.3962 

0 

0 

0.75 

SFIOO 

0.795 

0,2509 

0 

0 

1 

SFlOl 

0.795 

0.2509 

0 

0 

0.75 

SF102 

0.795 

0.5415 

0 

0 

1 

SF103 

0.795 

0.5415 

0 

0 

0.75 

SF104 

0.795 

0.3962 

4 

0 

1 

SF105 

0.795 

0.3962 

4 

0 

0.75 

SF106 

0.795 

0.2509 

4 

0 

1 

SF107 

0.795 

0.2509 

4 

0 

0.75 

SF108 

0.795 

0.5415 

4 

0 

1 

SF109 

0.795 

0.5415 

4 

0 

0.75 

SFllO 

0.795 

0.3962 

8 

0 

1 

SFlll 

0.795 

0.3962 

8 

0 

0.75 

SF112 

0.795 

0.2509 

8 

0 

1 

SF113 

0.795 

0.2509 

8 

0 

0.75 

SF114 

0.795 

0.5415 

8 

0 

1 

SF115 

0.795 

0.5415 

8 

0 

0.75 

SF116 

0.795 

0.3962 

12 

0 

1 

SF117 

0.795 

0.3962 

12 

0 

0.75 

SF118 

0.795 

0.2509 

12 

0 

1 

SF119 

0.795 

0.2509 

12 

0 

0,75 

SF120 

0.795 

0.5415 

12 

0 

1 

SF121 

0.795 

0.5415 

12 

0 

0.75 

SF122 

0.795 

0.3962 

0 

4 

1 

SF123 

0.795 

0.3962 

0 

4 

0.75 

SF124 

0.795 

0.2509 

0 

4 

1 

SF125 

0.795 

0.2509 

0 

4 

0.75 

SF126 

0.795 

0.5415 

0 

4 

1 

SF127 

0.795 

0.5415 

0 

4 

0.75 

SF128 

0.795 

0.3962 

0 

8 

1 

SF129 

0.795 

0.3962 

0 

8 

0.75 

SF130 

0.795 

0.2509 

0 

8 

1 

SF131 

0.795 

0.2509 

0 

8 

0.75 

SF132 

0.795 

0.5415 

0 

8 

1 

SF133 

0.795 

0.5415 

0 

8 

0.75 

SF134 

0.795 

0.3962 

0 

12 

1 

SF135 

0.795 

0.3962 

0 

12 

0.75 

SF136 

0.795 

0.2509 

0 

12 

1 

SF137 

0.795 

0.2509 

0 

12 

0.75 

SF138 

0.795 

0.5415 

0 

12 

1 

SF139 

0.795 

0.5415 

0 

12 

0.75 


100 






/ 




Appendix M. Results of Model Scale Experiments Performed 
at MIT to Determine the Restoring Forces and Moments due to 


Body Angle 

Shaded cells represent values greater than 1.15 standard deviations away from the mean 


Test# 

Subm. 

Spd 

Pitch 

Yatv 

X 

Y 

Z 

L__ 

Roll 

Pitch 

Yaw 

Normalized 
Sway Force 

Normalized 
Heave Force 

Normalized 
Yaw Moment 

Normalized 
Pitch Moment 

SF101 

252 

75 

0 

0 



gtgrgiii 

gliTigBH 

^■EESBi 

-0.0368 

0 

0 


0 

SF107 

252 

75 

0 

4 

KlidMcl 


tftlskUd 



0.0647 



EKEEB 


SF113 

252 

75 

0 

8 





-0.033 

0.1305 

EBEE^ 




SF119 

252 

75 

0 

12 


mW:Vkk\ 

giKtEH 


0.0022 


1.799 

BEHH 


BBBBH 

SF125 

252 

75 

4 

0 





0.0337 

0.0086 


BEE^ 


■■eezb 

SF143 

252 

75 

4 

0 



rntWiftm 

■EESS 

0.013 

0.0121 


EKES^ 


BBEESS 

SF131 

252 

75 

8 

0 



g*»i;lEl 


■■■CE^ 

0.0145 


EBEEiS] 


IBBESIi^EI 

SF137 

252 

75 

12 

0 



BUMld 

gplilcggl 

0.1141 

0.0176 




0.1557 

SF99 

398 

75 

0 

0 



giltfegid 

WtItUdd 


gBEESSQ 

HHHIKl 

0 


^^EBBuj 

SF105 

398 

75 

0 

4 



gigEwa 

tmiTx.ru\ 

HECEE)^ 

0.0572 

EEESSl 


■BEEEB 


SF111 

398 

75 

0 

8 



giT;ETgl 

MilikliVJ 

-0.044 


0.9341 




SF117 

398 

75 

0 

12 




glililiUd 

-0.0279 

EEXHEl 

1.5408 


BBEBES 


SF123 

398 

75 

4 

0 


BilEUn 

gig»EWH 


-0.003 

-0.0077 


-0.184 


BHEEIB 

SF129 

398 

75 

8 

0 

ggigni 


gigbEE! 

HESS] 

0.0613 

IEEEiES 


-0.719 


EBSES 

SF135 

398 

75 

12 

0 

■»nw«ia 


gi^»I;ll 


0.0867 

-0.0135 


-0.9307 


EHESEB 

SF103 

543 

75 

0 

0 

■EiEQ 




-0.0566 

EEEB 

0 

0 


0 

SF109 

543 

75 

0 

4 




M»MgH 

EHEE^ 

^EEE^ 

0.4138 


^^ee^ 


SF115 

543 

75 

0 

8 





-0.0898 

0.1038 

i^ESS 


0.0973 


SF141 

543 

75 

0 

8 





EHEE^ 

^I^HQSEQ 

EKES 


0.1275 


SF121 

543 

75 

0 

12 

1.4662 

0.8562 

0.0813 

0.0017 

-0.0339 

0.1749 

■miES 

HHBBi 

0.1684 

^^BBH 

SF127 

543 

75 

4 

0 





0.0111 

EEEB 


ibeeh 


ehe^ 

SIEKH 

543 

75 

8 

0 



gtlMkl 


HHKE^E] 

-0.0012 


ibeeeb 


f^B^EXES 


543 

75 

12 

0 


gtHElHU 

g*mdEi 

ItliMld 

■EKES] 

-0.004 


■EEE 


bbexe^ 


252 

■■1S3 

0 

0 

K;Igc!;Sl 

g*lMM.i 

g*r$Ei 



-0.0408 


0 

BBBIH] 

0 

SF106 

252 

IHIiSl 

0 

4 



guMbii 

gfililiEH 

-0.0736 

0.0959 

HES^ 


0.1367 


SF112 

252 


0 

8 





-0.0857 

0.199 

■EEQS 


BE^^ 


SF118 

252 

100 

0 

12 


■ntnsia 



-0.0386 

■ESBSj 

eibes 


BBEES30 


SF124 

252 

100 

4 

0 




HESSU 

0.0335 

0.0038 


-0.648 


BEEE^ 

SF130 

252 

100 

8 

0 





0.0674 

0.0062 


BHEEu 


BEEXI^ 

SF142 

252 


8 

0 





0.0661 

0.0052 


BHEEiE3 


BBEXB& 

SF136 

252 

100 

12 

0 

1MHI 

giWg<Zi 


■EE^ 

0.1288 

HEEESEB 


^ebes 


EBEXESS 

SF98 

398 

100 

0 

0 

■ElcIWH 




-0.1263 

-0.0457 

HKHEj 

0 


BBI^BD 

SF104 

398 

100 

0 

4 


gtKitwa 

gitZIiH 

MtXtkkbl 

-0.1217 


lEESESS 


0.1116 



■■ES! 

100 

0 

8 

mKMA 


gtgggfn 

gtXiEIH 

fmUQQQQ 

HBKKE&l 

^EESS 


EK^D 





0: 

12 



■fiiad 

piilitwm 

-0.0714 

0.2824 

2.2185 

^^BBB 

BBEE^ 

^^BHE 


398 


4! 

0 



giltViiVi 


-0.0199 

-0.0161 


HBEESi 


BBO&l 


398 

100 

8j 

0 


gtKltUM 

g«ig:m 

gfigiliEH 


-0.0037 


-1.2357 


BBEXESI 


398 

100 

12 

0 


gigbgM 



0.1172 

EEES 


-1.5815 


BBCEES 

SF102 

543 

100 

0: 

0 




Mili'i™ 

-0.1139 

0 

0 

0 

0 

0 

SF108 

543 

100 

0 

4 

13.1946 

0.6244 

-0.1566 

0.0198 

-0.0778 

0.0915 

EKEni 


0.0915 


SF114 

543 

100 

0 

8i 


gigrgigi 



■■■EliEE] 

0.1744 

EKB3Q 


0.1744 


SF140 

543 

100 

0 

8| 


gtl;g!kgl 



■■BEUD 






SF120 

543 

■■IE13 

0 

12 


wsKym 



EESBiB 

0.2688 

■BBsE 




SF126 

543 


4 

0 

1.5582 

-0.3684 

-0.1921 

0.0073 

0.0015 

. -0.0231 


BBEEiSS 


■EOXEil 


543 

100 

8 

0 

■nCTHiH 




0.0606 

-0.0079 


-0.8489 


0.1745 


543 

100 

12 

0 





0.1005 

-0.0276 


-1.0989 


0.2144 


€ 


Units for Interpreting the Results of Model Scale Experiments Performed at MIT to Determine the 
Restoring Forces and Moments due to Body Angle 


Quantity 

Units 

Depth 

cm 

Submergence 

mm 

Speed 

Cm/s 

Angle 

Degrees 

Forces 

Newtons 

Moments 

N-m 


101 


























































